SMOOTH TYPE II BLOW UP SOLUTIONS TO THE FOUR DIMENSIONAL 
ENERGY CRITICAL WAVE EQUATION 
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Abstract. We exhibit C°° type II blow up solutions to the focusing energy critical wave equation 
in dimension N = 4. These solutions admit near blow up time a decomposiiton 

u(t,x)= (Q + e (t))(-r^) with ||e(t),&e(t)||ai x£a «l 

where Q is the extremizing profile of the Sobolev embedding H 1 —¥ L 2 , and a blow up speed 

A(t) = (T — f) e -VH°^-*)l(i+°(i)) as t T. 



1. Introduction 

1.1. Setting of the problem. We deal in this paper with the energy critical focusing wave equation 



I d tt u - An - f(u) = with f(t) = t^-a, a - 

\ (« J 5tt*)|t=o = (wo,wi) J GMxM^. 

in dimension 

N = 4. 

This is a special case of the nonlinear wave equation 

d a u-Au-f(u) = (1.2) 
which since the pioneering works by Jorgens [10] has attracted a considerable amount of works. 

JV+2 

For the energy critical nonlinearity f(u) = ±t N ~ 2 , the Cauchy problem is locally well posed in 
the energy space H 1 x L 2 and the solution propagates regularity, see for example Sogge [35] and 
references therein. Recall that in this case, (1.2) admits a conserved energy 

If n If n N ~ 2 f 2N 

E(u(t)) = E(u , ui) = - J {d t uf + -J |Vn| 2 T J 
which is left invariant by the scaling symmetry of the flow: 

A— A X 2 

Global existence in the defocusing case was proved by Struwe |38| for radial data and Grillakis [U] 
for general data. For focusing nonlinearities, a sharp threshold criterion of global existence and 
scattering or finite time blow up is obtained by Kenig and Merle |14j based on the solitonic solution 
to ( O ): 

JV-2 

Q{r) = ( - ^5— 1 2 (1-3) 

V 1 + N(N-2) J 

which is the extremizing profile of the Sobolev embedding H 1 — > L 2 * . Indeed, for initial data (no, u\) 
such that E(uq,ui) < E(Q,0), those with || Vuq || z, 2 < II^QIIl 2 have global solutions and scatter, 
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while those with [|Vito||£2 > ||VQ||l 2 l ea d to finite time blow up. 

Note that like in the works by Levine [20], see also Strauss [37], and as is standard in a nonlinear 
dispersive setting, blow up is derived through obstructive convexity arguments, see also Karageorgis 
and Strauss [IT] for refined statements near the soliton Q. However, this approach gives very little 
insight into the description of the blow up mechanism and the description of the flow even just near 
the ground state soliton Q is still only at its beginning. 



1.2. On the energy critical wave map problem. There is an important litterature devoted to 
the construction of blow up solutions for nonlinear wave equations, see e.g. Alinhac [JJ, and Merle 
and Zaag [28J , [29J for the study of the ODE type of blow up for subcritical nonlinearities. For energy 



d tt u-d rr u r - — = 0, (1.4) 



critical problems like (1.1), recent important progress has been made through the study of the two 
dimensional energy critical corotational wave map to the 2-sphere: 

d r u k 2 sin 2u 
272 

where k £ N* is the homotopy number. The ground state is given there by 

Q{r) = 2tan" 1 (r fc ). 

After the pioneering works by Christodoulou, Tahvildar-Zadeh [4], Shatah and Tahvildar-Zadeh [36] 
and Struwe [39] and their detailed study of the concentration of energy scenario, the first explicit 
description of singularity formation for the k = 1 case is derived by Krieger, Schlag and Tataru |17j 
who construct finite energy finite time blow up solutions of the form 

u(t,x) = (Q + e)(t,^r) with ||e(t),5te(*)llHi x ^< 1 ( L5 ) 
with a blow up speed given by 

\(t) = (T-ty, vi/>~, 

see also |19j . The spectacular feature of this result is to exhibit arbitrarily slow blow up regimes 
further and further from self similarity which would correspond to the -forbidden, see [39J- self 
similar law 

A(t)~r-t. (i.6) 

Numerics suggest [3] that this blow up scenario is non generic and corresponds to finite codimensional 
manifolds. After the pioneering works [34j for large homotopy number k > 4, Raphael and Rodnianski 
|31] give a complete description of a stable blow up dynamics which originates from smooth data 
and for all homotopy number k > 1. The blow up speed obeys in this regime a universal law which 
depends in an essential way on the rate of convergence of the ground state Q to its asymptotic value 

ir — Q ~ — as r — > oo, 



and indeed the stable blow up regime corresponds to a decomposition (1.5) with the blow up speed 



c k T ~ l t for k>2, 

\(t)~{ |io g( r-i) p-2 n 7 ) 

{T _ t)e -VW^y\ for k = l. 

Note that this work draws an important analogy with another critical problem, the L 2 critical 
nonlinear Schrodinger equation, where a similar universality of the stable singularity formation near 
the ground state is proved by Merle and Raphael in the series of papers [23], [21], [30], [2S], [26] . 
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1.3. Statement of the result. For the power nonlinearity energy critical problem ( |1,1[ ), there has 
been recent progress towards the understanding of the flow near the solitary wave Q. In [15], Krieger 
and Schlag construct in dimension N = 3 a codimension one manifold of initial data near Q which 
yield global solutions asymptotically converging to the soliton manifold. The strategy developed 
by Krieger, Schlag, Tataru in [T7J for the wave map problem has been adapted in |18| to show in 
dimension N = 3 the existence of finite energy finite time blow up solutions of the form 

1 x 

u(t,x) = — ^ — (Q + e)(t,—rpr) with \\e(t), Qte(t)||^i xL a < 1 
A 2 (t) 

and with a blow up speed given by 

\(t) = (T-ty, Vz/>| (1.8) 

The quantization of the energy at blow up for small type II blow up solutions in dimension N 6 {3, 5} 
is proved in [6], [7J in the radial and non radial cases. In particular, for radial data, if T < +oo and 

sup [\Vu{t)\\ 2 + d t u\ 2 L2 ] < \VQ\ 2 L2 +a*, a* < 1, 
te[o,T] 

then there exists a dilation parameter A(t) — > as t — > T and asympotic profiles (u*,v*) 6 H l x I? 
such that 

\u(t,x) ^ — Q(^r),d t u(t) ) -> (u*,t>*) in H 1 x L 2 as t -> T, 

see [27] for related classification results for the L 2 critical (NLS). 

These works however leave open the question of the existence of smooth type II blow up solutions. 
We claim that such smooth type II blow up solutions can be constructed in dimension N = 4 as the 
formal analogue of the singular dynamics exhibited by Raphael and Rodnianski |31j for the wave 
map problem in the least homotopy number class k = 1. The following theorem is the main result 
of this paper: 

Theorem 1.1 (Existence of smooth type II blow up solutions in dimension N = 4). Let N = 4. 

Then for all a* > 0, there exist C°° initial data (uq,u\) with 

E(u ,ui) < E(Q,0)+a* 



such that the corresponding solution to the energy critical focusing wave equation (1.1) blows up in 
finite time T = T(uo,Ui) < +oo in a type II regime according to the following dynamics: there exist 
(u*,v*) S H 1 x L 2 such that 

(u(t,x) ^ — Q(-^),d t u(t) ) ->• (u*,v*) in H 1 x L 2 as t -»■ T (1.9) 

V A-r(f) A W / 

with a blow up speed given by 

X(t) = (T — i) e -Vl lo sCr-*)|(i+°(i)) as t ^ T _ ( L10) 

Comments on the result 

1. On the smoothness of the initial data: An important feature of Theorem |1.1| is to exhibit a 
new blow up speed which is valid for C°° solutions. Indeed, while the Krieger, Schlag, Tataru |18j 
approach provides a continuum of blow up speeds, the exact regularity of the obtained solutions is 
not known, which is an unpleasant consequence of their construction scheme. In fact, it is expected 
that C°° initial data should lead to quantize blow up rates hence breaking the continuum of blow up 



speeds (1.8), we refer to [2] for a related discussion in the context of the energy critical harmonic 
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heat flow. Hence we expect the blow up rate (1.10) to correspond to the minimal type II blow up 



speed of smooth solutions with small super critical energy. Such a general lower bound on blow up 
rate in the spirit of the one obtained by Merle and Raphael for the I? critical NLS [30], [26J is an 
open problem. The construction of excited blow up solutions with other speeds and C°° regularity 
also remains to be done. This problematic is related to the understanding of the structure of the 
flow near Q which is still at its beginning. 



2. On the codimension one manifold: The proof of Theorem 1.1 involves a detailed description of 
the set of initial data leading to the type II blow up with speed (1.10). Indeed, given a small enough 
parameter bo > and a suitable deformation Qi, of the soliton with 

Qb Q as b ->• 

in some strong sense, we show that for any smooth and radially symmetric excess of energy 

\\V0,m\\H2xHl $ 



log(6o)| 



we can find d_|_(6o, 770, 771) € K such that the solution to (1.1) with initial data 



uo = Qb +Vo + d + tp, ui 



N 



Qbo + y • ^Qb + m 



blows up in finite time in the regime described by Theorem |1.1[ Here ip is the bound state of 
the linearized operator close to Q and generates the unstable mode, we refer to Definition 3.4 and 



Proposition 3.5 for precise statements. Hence the set of blow up solutions we construct live on a 



codimension one manifold in the radial class in some weak sense. Following [TS], [IB], the proof that 
this set is indeed a codimension one manifold relies on proving some Lipschitz regularity of the map 
(bo,i]OiVi) ~^ d+(bo,i]o,r]i), and in particular some local uniqueness to begin with. The analysis in 
[16] shows that this may be a delicate step in some cases. Our solution is constructed using a soft 
continuous topological argument of Brouwer type coupled with suitable monotonicity properties in 
the spirit of Cote, Marte and Merle [5], and in other related settings, see e.g. Martel [21], Raphael 
and Szeftel [32], this strategy has proved to be quite powerful to eventually achieve strong unique- 
ness results. This interesting question in our setting will require additional efforts and needs to be 
adressed separately in details. 



3. Extension to higher dimensions: We focus onto the case of dimension N = 4 for the sake of 
simplicity, and our main objective is to provide a robust framework to construct C°° type II blow up 



solutions. However, following the heurisitic developed in [31], the blow up speed (1.10) corresponds 
to the k = 1 case in (1.7), and we similarly conjecture in dimension N > 5 the existence type II 



finite type blow up solutions close to Q with blow up speed 

T — t 

A(i) ~ c N - 



|log(T-t) 



1 

JV-4 



Note from (1.3) that the higher the dimension, the fastest the decay of the ground state Q, and 



this should avoid some difficulties which occur only in low dimension like in [31] for large homotopy 
number k > 4. We expect the strategy developed in this paper to carry over to the case N = 5,6, 
but the extension to large dimension will be confronted in particular to the difficulty of the lack 
of smoothness of the nonlinearity. Let us also insist onto the fact that the case N = 4 is in many 
ways the more delicate one in terms of the strong coupling of the main part of the solution and 
the outgoing tail due to the slow decay of Q, which results in the somewhat pathological blow up 
speed (1.10). This comment becomes even more dramatic in dimension N = 3 where we expect 



our analyis to be applicable to construct C°° type II blow up solutions, but this seems to require a 
slightly different approach. 

1.4. Strategy of the proof. Let us briefly summarize the strategy of the proof of Theorem |1.1[ 
step 1 Approximate self similar solution. 



Let -D, A denote the differential operators (1.18). Exact self similar solutions to (1.1) of the form 

u(t, x) = -j^- 2 — Q b [ ] with b = -A t 
A~a~(i) \ A W/ 

where Q b satisfies the self similar equation 

AQ b -b 2 DAQ b + Q 3 b = (1.11) 

are known to develop a singularity on the light cone y = = g leading to an unbounded Dirichlet 
energy ||VQb||^2 = +oo, see Kavian, Weissler [12j. We therefore assume < b >C 1 and consider a 
one term expansion approximation 

Qb = Q + b 2 T ± 



which injected into (1.11) yields at the order b 2 : 

HT\ = -DAQ. (1.12) 

Here H is the linearized operator close to Q given by 

N + 2 4 

H = -A- ] ^Q^2. (1.13) 

The spectral structure of H is well knwon in connection to the fact that Q is an extremizer of the 
Sobolev embedding H 1 — > L 2 , and in the radial sector, H admits one non positive eigenvalue with 
well localized eigenvector ip: 

H^ = -C^P, C>0, (1.14) 

and a resonance at the boundary of the continuum spectrum generated by the scaling invariance of 
(0): 

H(AQ) = 0, AQ(r) ~ as r -> +oo. (1.15) 

In order to solve (1.12), we first remove the leading order growth in the exact solution T\ = \ \y\ 2 Q 
which is consequence of the flux computation: 



{DkQ,kQ) = \ lim y 4 |AQ| 2 >0 (1.16) 
Z y— >+oo 



due to the slow decay of Q in dimension iV = 4 from (1.3). For this, we solve 

HTi = -DAQ + c&AQl <i with c b = ~ ^Ti as b ~> °- 

The purpose of this construction is to yield after a suitable localization process an o(b 2 ) approximate 
solution to the self similar equation ( |1.11[ ) which dominant term near and past the light cone is still 
given by Q itself in the sense that: 

6 2 |Ti| <Q for y > i 
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This identifies Q as the leading order radiation terrrj^J 
step 2 Bootstrap estimates. 



We now roughly consider initial data of the form 

uo = Qb + d+ip + 770, «i = M-Q&o + 7/i, with \d+\ + ||t/o, r?i||#2 x #i < b%, 
and introduce a modulated decomposition of the flow 



(1.17) 



u(t, x) 



1 



iV-2 , 



A(i) 



6(f) 



A*. 



Here we face the major difference between the power nonlinearity wave equation ( 1.1 ) and the critical 



wave map problem (1.4) which is the presense of a negative eigenvalue in the first case (1.14) for the 



linearized operator H close to Q. This induces an instability in the modulation equations for b, A 
which is absent in the wave map case, leading to stable blow up dynamics. However, we claim that 



the ODE type instability generated by (1.14) is the only instability mechanism 



The situation is conceptually similar to the one studied in [5] where multisolitary wave solutions 
are constructed in the supercritical regime despite the presence of exponentially growing modes for 
the linearized operator which are absent in the subcritical regime. We adapt a similar scheme of 
proof which does not rely on a fixed point argument to solve the problem from infinity in tim^J but 
by directly following the flow for any initial data of the form (1.17). This reduces the full problem 



to a one dimensional dynamical system for which a classical clever continuity argument yields the 
existence of d+(&o> ?/o, Vi) such that the unstable mode is extinct, see section [5j 
The key is hence to control the flow under the a priori control of the unstable mode, and here 
we adapt the technology developed in [31] which relies on monotonicity properties of the linearized 
Hamiltonian at the H 2 level of regularity. However, the analysis in [31] heavily relies on the existence 
of a decomposition of the Hamiltonian 

H = A* A, A = -d y + V{y) 

which is central in the proof of the main monotonicity property and is lost in our setting. This forces 
us to revisit the approach in several ways, and to rely in particular on fine algebraic properties of 
the flovJ^ near Q and coercitivity properties of suitable quadratic forms in the spirit of [22] . [23], see 
Lemma 4.7 which remarkably turn out to be almost explicit thanks to the formula (1.3). We are 
eventually able to find d-f(&o, rjo, r/i) for which to leading order 



-c b b 2 ~ 



ds 
dt 



1 

A' 



\d+ \ + 



\dyy £ \ 



L 2 



<C b 2 



2|log6|' 

which reintegration in time yields finite time blow up in the regime described by Theorem |1.1[ 



1.5. Notations. We define the differential operators: 
N-2 



A/ 



N 



f + y-Vf (IT 1 scaling), Df=-f + yVf (L 2 scaling). 



(1.18) 



Denoting 



/r+oo 
f9 = J Q f(r)g(r) 



* N - l dr 



^see [31] for a further discussion on this issue and the role played by the non vanishing Pohozaev integration ( 1.16 1 
2 after renormaliza tion o f the ti me 
3 see in particular k.23|, k.38 1 
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the L 2 (R N ) radial inner product, we observe the integration by parts formula: 

(Df,g) = -(f,Dg), (Af, g) + (Ag, f) = -2(/, g). (1.19) 
Given / and A > 0, we shall denote: 

and the space rescaled variable will always be denoted by 

r 

y= x- 

We let x be a smooth positive radial cut off function x( r ) = 1 for r < 1 and %(r) = for r > 2. For 
a given parameter B > 0, we let 

XB(r) = X (£) • (1-20) 

Given 6 > 0, we set 

*4 a = !M. (1 . 21) 

To clarify the exposition we use the notation a <b when there exists a constant C with no relevant 
dependency on (a, b) such that a < Cb. In particular, we do not allow constants C to depend on the 
parameter M except in Appendix A. 

Aknowledegments: The authors would like to thank Igor Rodnianski for stimulating discussions 
about this work. P.R is supported by the French ANR Jeune Chercheur SWAP. 



2. Computation of the modified self-similar profile 

This section is devoted to the construction of an approximate self-similar solution Qt which de- 
scribes the dominant part of the blow up profile inside the backward light cone from the singular 
point (0, T) and displays a slow decay at infinity which is eventually responsible for the modifications 
to the blow up speed with respect to the self similar law. The key to this construction is the fact 
that the structure of the linearized operator H close to Q is completely explicit in the radial sector 
thanks to the explicit formulas at hand for the elements of the kernel. 

We introduce the direction 

$ = DAQ (2.1) 

which displays the cancellation 

My)\ Z (2-2) 

and the crucial nondegeneracy which follows from the Pohozaev integration by parts formula: 

(d>,AQ)= Urn (ly 4 \AQ\ 2 ) = 32 > 0. (2.3) 

Proposition 2.1 (Approximate self-similar solution). Let M denote a large enough constant. Then 
there exists b*(M) > small enough such that for all < b < b*(M), there exists a smooth radially 
symmetric profile T\ satisfying the orthogonality condition 

(Ti, X m$) = (2.4) 

such that 

Pb 1 =Q + Xb 1 1?T 1 (2.5) 
is an approximate self similar solution in the following sense. Let 

* Bl = -AP Bl + b 2 DAP Bl - f(P Bl ), (2.6) 
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then for all k > 0, < y < ^. 



d k Tx 
dy 

d k dP Bl 



dy k db 



~ l + y k 



< bly<2B 1 



1 log(M) + |log(l+y)| 

log&| ^<v<^ b 2 y 2 \\ogb\ v>^ 1 + y 2 



l + \log(by)\ 1 
|log6| 

1 + |log(6y)| 



1+y* 



Sq + 



1 



log&| 2<y<^ 6 2 y 2 |log | y>^ 



+ 



log(M) + |log(l + y)| 
1 + y 2 



and, for all k > 0, y > 0, 
dy 

6 4 ri + |log(6y)| 



< 

~ 1 + y fc 



Sq + 



logb| 2<y<Jl 6 2 y 2 |log6| 2B l>?/> 



Sfl + 



log(M) + |log(l + 
1 + y 2 



Lj/<2Bi 



+ (1 _|_ y4+fc) 1 J/>- B i/ 2 ' 

/or some constant 



Cb 



2|log6| 



1 + 



|log6| 



Proof of Proposition |2.1 
Step 1 Inversion of ii\ 

The first green function of -ff is given from scaling invariance by 



2 1 + 



AQ(y) 



which admits the following asymptotics: 

o*(AQ) 



AT-2 



2/ 

N(N-2) 



N(N - 2) 



V k > 0, 



dy* 



(y) 



O(l) as y^0, 

G(y -(iV-2+fc)) 



as y — )• oo 



Let now 



r(y) = -AQ(y) 



ds 



-JV-1 



(Ag) 2 ( s )' 



(2.7) 



2i 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



be another (singular at the origiiQ element of the kernel of H which can be found from the Wronskian 
relation: 

T'AQ - r(AQ)' ~ ' 



JV-1 ■ 



From this we easily find the asymptotics of r( fc ) for any integer k: 



d k T 
dy 



k (y) 



0(y~ k ) as y — > oo. 



A smooth solution to Hw = F is given by: 

w(y) = T(y) f V F(s)AQ(s)s N ~ 1 ds - AQ(y) [ F{s)r{s)s N ' 1 ds. 
Jo Jo 



(2.13) 



(2.14) 



Note that F must be smooth at y — \J N(N — 2) where AQ vanishes from the radial ODE HV — 



We now look for a solution to the self similar equation in the form Q + b 2 T\. This yields: 



: -AQ b + b 2 DAQ b - f(Q b ) 
= b 2 {HT x + DAQ) + b 4 DATi - [f(Q + 6 2 Ti) - f(Q) - b 2 f{Q)T x } . 
Step 2 Computation of T\. 



(2.15) 



Thanks to the anomalous decay (2.2), we chose T\ solution to 

( HT\ = F = -DAQ + c fe xso AQ, 



with Cf, chosen such that: 



(F,AQ) = 



(2.16) 
(2.17) 



i.e. from Pohozaev integration by parts formula, see (1.21) and (2.3) , 

Cb 



(DAQ,AQ) __ l lim^ +00 /|AQ(j/)| 2 
( X soAQ,AQ) 2 Jx^olAQI 2 



1 



2|log6| 



1 + 



1 



|log&| 



as 6—t-O. 



This yields (2.10). Following (2.14), we first consider 



f 1 (y)=T(y) F(s)AQ(s)s 3 ds - AQ(y) I F(s)T(s)s 3 ds 



(2.18) 



The smoothness of T\ at the origin follows from (2.18) together with elliptic regularity from (2.16). 
We now examine the behavior of T% at large y. 



We first observe that, from the orthogonality (|2.17|): 

Ti(y) = 



r(y) 



+oo 



F{s)AQ(s)s 3 ds + AQ(y) / F(s)T(s)s 3 ds 



Hence, from the degeneracy |DAQ| = 0(y 4 ), this yields that, for ^ < y < -jk: 



\Ti(y)\ 



< 



(1 + S 4 )(1 + S 2 ) 



ds + 



l + s 4 



ds + \cb\ 



Bo s 3 



1 + S 2 



zds 



< Mi+J 



1+2/ 2 



6 2 y 2 |log6| ' 



(2.19) 



similarly, for 1 < y < 

\Ti(y)\ ~- 



< 



+ 



+oo 

1 



+oo 



F(s)AQ(s)s 3 ds + AQ(y) / F(s)T(s)s 3 ds 



(l + s 4 )(l + s 2 



-ds + |c 6 | 



So 



2/ s 3 



1 + S' 



;ds + \c b \ 



y s 3 



l + s 2 



(l + s 2 ) 2 
ds 



ds 



< l + |log(6y)| | |log(l + y)| 



|log6| 



i + y 2 



(2.20) 
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We now choose thanks to (2.3): 



T 1 (y) = f 1 {y)-cAQ with c 



Cfi,XAf$) 
(XM$, AQ) 



so that the orthogonality condition (I2.4I) is fulfilled. We note that so that the bounds ( 2.19 ) and ( 2.20 ) 



ensure that c remains bounded by log(M) uniformly in M and b, provided b is chosen sufficiently 
small w.r.t. M. 



This yields (2.7) for k = 0, the other cases follow similarly. 



Step 4 Estimate on ^ Bl and d b ^ Bl 



We now cut off the slow decaying tail T\ according to (2.5) and estimate the corresponding error 



to self similarity ^ Bl given by (2.6). 
We compute: 



tffl! = b 2 XBi {HTi + DAQ) 

+ b 2 [-2 XBl T[ - TxA XBl + (1 - XbJDAQ + b 2 DA{xB,T x )] 

- [f{Q + b\ Bl T{)- f(Q)-XB 1 f'(Q)T 1 \. 

Outside the support of xb 1 we have thus ^ Bl = b 2 DAQ. On the other hand, in dimension N = 4, 
we have the Taylor expansion : 

HQ + b 2 XBl Ti) - HQ) - XbJ{Q)Ti = b\ 2 Bl T 2 {y) f\l - r)(Q(y) + rb^T^dr. 

Jo 



We thus estimate from (2.7), (2.15), (2.16) and the degeneracy (2.2) for y < 2B\ : 

1 



<b'l 



+ 



Ti 



+ 



y>B 1 / 2 y i + y i 1+y2 . 1 + y4 



+ ft^DA^Ti)! + 6 4 |2f / (1 - r)|Q(y) + T&M^dT. 

JO 



(2.7) now yields (2.9) for fc = 0. Further derivatives are estimated similarly thanks to the smoothness 



of the nonlinearity. We emphasize here that, given B > large, we have 1/(1 + y) < 1/B < 1/(1 + y) 
on the support of x' B , so that differentiating xb acts as a multiplication by 1/(1 + y). Furthermore, 
there holds 1/ B\ = o{b) so that we can always dominate 1/(1 + y) by b on the support of x' Bl - 
Finally, we compute dbPB 1 from (2.5). 

To this end, we note that d b c b = 0(l/6|log(6)| 2 ) when b — > so that the source term for T\ in 
(2.16) satisfies 



d b F 



AQ 



where p{z) = zx'(z) £ C£°(0,oo) and we keep the convention for function dilation. Hence, the same 
arguments as for T\ enable to show that d b T\ and then d b T\ satisfy the estimates: 



dy k 



(y) 



< 



1 



6(i + v k ) 



l + |log(6y)| l_ 

2<y<^- + b 2 y 2 \\ogb\ ~v> 



|log6| 



1 ^ + 

2 



l + |log(l + y)| 



i + y 2 



Finally, we compute from (2.5) 



d b P Bl = 2b XBl T 1 + 6 2 ^log(Bi)p Bl Ti + b 2 XBl d b Ti. 



(2.21) 
(2.22) 



This decomposition together with (2.7) and the previous computation yield (2.8) 



This concludes the proof of Proposition 2.1 
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3. Description of the trapped regime 

We display in this section the regime which leads to the blow up dynamics described by Theorem 

o 



3.1. Modulation of solutions to (1.1). Let us start with describing the set of solutions among 



which the finite time blow up scenario described by Theorem |1.1| is likely to arise. We recall from 



(1.14) that ip denotes the bound state of H with eigenvalue — £ < 0. The following lemma is a 
standard consequence of the implicit function theorem and the smoothness of the flow, see Appendix 

m 

Lemma 3.1 (Modulation theory). Let M be a sufficiently large constant to be chosen later and 
< 6o < Oq{M) small enough. Let (770, r/i, c£+) satisfying the smallness condition: 

b 2 

\d+\ + ||r ? o I V7?o,r?i + 6o(l-XB 1 (6 ))AQ,Vr ? i||^ lx ^ lxi . 2xi . 2 < ^J^, (3.1) 
then, there exists a time Tq such that the unique solution u e C 2 ([0, To]; L 2 (R N )) nC([0, To]; H 2 (Mr)) 



to (1.1) with initial data : 

u o = PBtlbo) +Vo + d+ip, ui = b AP Bl{bQ) + 771, (3.2) 
admits on [0, To] a unique decomposition 

u(t) = (P Bl{m + e(t)) x{t) (3.3) 

with 

LAG C 2 ([0,To],R;) with 

Vie [0, T ], (e(i),XAf*) = and b(t) = -X t ; (3.4) 
2. there holds the smallness: 

\\Ve(t)\\ L 2 < 6 |log6 | \b(t) - b \ + \X(t) - 1| + \\V 2 e(t)\\ L , < V t e [0,T ]. (3.5) 

|logo | 

Remark 3.2. Recall that the slow decay of Q and the choice of P Bl induces an unbounded tail of 
Tb 1 in the energy norm, and more specifically [|AQ[|^ = +oo, hence the need for the compensation 



in the norm for the time derivative in (3.1). 



3.2. Decomposition of the flow and modulation equations. Considering initial data satisfying 



the assumption of the above lemma, we now write the evolution equation induced by (1.1) in terms 



of the decomposition (3.3). Let 

1 / „ s / r 



where b = —Xt- Let us derive the equations for w and e. Let 



°®=Jlw) (3 ' 7) 

be the rescaled time. We shall make an intensive use of the following rescaling formulas: 

u(t,r) = J ^ =I v(s,y), y= r -, | = i (3.8) 

8 t u = \ (d s v + bAv) x , (3.9) 
A 

d tt u = [d 2 v + b(d s v + 2Ad s v) + b 2 DAv + b s Av] x . (3.10) 
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In particular, we derive from (1.1) the equation for e: 

8 2 e + H Bl e = -tt Bl - b s AP Bl - b(d s P Bl + 2kd s P Bl ) - d 2 s P Bl 
- b{d s e + 2Ad s e) - b s Ae + N(e) 



(3.11) 



where, implicitly, B\ = B\{b{t)) and H Bl is the linear operator associated to the profile P Bl 

H Bl e = -Ae + b 2 DAe-f'{P Bl )e, (3.12) 

and the nonlinearity: 

N(s) = f(P Bl + e) - f(P Bl ) - f(P Bl )e. (3.13) 
Alternatively, the equation for w takes the form: 

d 2 t w + H Bl w = - [d 2 {P Bl ) x - A(P Bl ) x - /((PbJa)] + N\(w) 

with 

H Bl w = -Aw- f'((P Bl ) x )w, (3.14) 
N x (w) = f((P Bl ) x +w)- f((P Bl )x) ~ f'((P Bl )x)w. (3.15) 



We then expand using (3.9), (3.10): 

d 2 (P Bl )x-A(p Bl ) x -f((p Bl ) x ) 

and rewrite the equation for w: 

1 



A 2 
1 



[d ss p Bl + b(d s p Bl + 2Ad s p Bl ) + b s AP Bl + v B ] 



^ [bAd s p Bl + b s AP Bl + i- B ] A + d t 



\{d s P Bl )x 



d 2 w + H Bl w = -^ [bAd s P Bl + b s AP Bl +* B ] X - d t 



j(d s P Bl )x 



+ N x (w). (3.16) 



For most of our arguments we prefer to view the linear operator H Bl acting on w in ( |3.16 ) as a 
perturbation of the linear operator H x associated to Q x . Then 

d?w + H x w = F Bl (3.17) 



with 



-2 [bAd s p Bl + b s AP Bl + * Bl ] A - d t 

-[f'(Q x )-f\(P Bl ) x )} w + N x (w) 



H x w = -Aw + f(Q x )w. 



(d s p Bl ; 



(3.18) 



3.3. The set of bootstrap estimates. At first, we fix some notations. We introduce the energy 
£(t) associated to the Hamiltonian H x : 



£(t) = X 2 J [(H x d t w, d t w) + (H x w) 2 ] . 



Given £ G (0,oo) the unstable eigenvalue, we set: 

1 



V- 



1 



and, we introduce the decomposition of the unstable direction 



{d s e, tp) 



a + (s)V + + d-(s)V- 



Let us denote: 



k+(s) = a + (s) + ^=(d b P Bl ,ip), K-(s) = a-(s) - ^=(d b P Bl , ?p) . 



(3.19) 

(3.20) 

(3.21) 
(3.22) 
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We note that the vectors V+, V- given by (3.20) yield an eingenbasis of 

1 

C 

and hence correspond respectively to the unstable and stable mode of the two dimensional dynamical 
system 

ds I C o ! 



which to first order in b is verified by the projection onto the unstable mode (e, i[>), see (4.57). The 
deformation term b s (df,PB 1 ,' t P) in ( 3.22 ) is present to handle some possible time oscillations induced 
by the 8 2 Pb 1 term in the RHS of (3.11) which cannot be estimated in absolute value but will be 
proved to be lower order. 

With these conventions, we may now paramaterize the set of initial data described by Lemma 



3.1 by a+ = k+(0), and then reformulate the initial smallness properties in terms of suitable initial 



bounds for e, see Appendix [A] for the proof which is standard. 

Lemma 3.3 (Inital parametrization of the unstable mode and initial bounds). Let M and bo be 

given as in Lemma 3.1 and denote by C{M) a sufficiently large constant. Then, given (770,771,04.) 
satisfying 

b 2 

l«+l + ll%,Vr / o,7 ?1 + 6o(l-XiJ 1 (b )) A Q' V7 ?illijixHixL2xL2 ^ \hgb ~\ , ( 3 ' 23 ^ 
there exists a unique d+ with \d+\ < 6g/|log(&o)| an d Tq > such that the unique decomposition 

U{t) = {PB t (b(t))+E)\(t) = (PBi(b(t)))\(t) +W{t), 
of the unique smooth solution u to (1.1) on [0, To] with initial data (3.2) satisfies the initialization 

K+ (0) = a+, (3.24) 
and the following smallness condition on [0, To] : 

• Smallness and positivity of b: 

< b{t) < 5b ; (3.25) 

• Pointwise bound on b s : 

• Smallness of the energy norm: 

\\(Vw(t),d t w(t) + M((i -XB l( 6(t)))AQ)) A(t) || L2xL 2 < Vb~o; (3.27) 

• Global H 2 bound: 

|f(t)|£C(M) W ; (3 ' 28) 



A priori bound on the stable mode: 

< (r,(M\\\ 

|log&(t)| 



«_(t)| < (C(M))* } b ®}* . (3.29) 



A priori bound of the unstable mode: 



|log&(i)| 

We may now describe the bootstrap regime as follows: 



K+(t)\ < 2 ,fy.. . (3.30) 
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Definition 3.4 (Exit time). Let K(M) denote some large enough constant. 



Given a+ € [— rfc^n > |i O g°b | L we Zei T(a+) 6e the life time of the solution to (1.1) with initial data 



(3.2), and Ti(a+) > be the supremum of T £ (0,T(a+)) such that for all t G [0, T] , t/ie following 
estimates hold: 

• Smallness and positivity of b: 

< 6(i) < 56 ; 

• Pointwise bound on b s : 



< K(M) 



|log6(t)| ; 



(3.31) 
(3.32) 



Smallness of the energy norm: 



bit) 



\\(Vw(t),d t w(t) + -^((1 - XB 1 {b(t)))&Q))Kt)\\v>xL» 



Global H 2 bound: 



A(i) 



\£{t)\<K{M)- 



< 



|log6(t)|2' 

• ^4 priori bound on the stable and unstable modes: 



M*)l <2 



|k_(<) | < (A"(M))b 



(3.33) 
(3.34) 
(3.35) 



|logb(t)|' ------- |log6(t)|' 

The existence of blow up solutions in the regime described by Theorem |1.1| now follows from the 
following: 



Proposition 3.5. There exists a+ G 



|log6o| ' |logbo| 
Ti(a+) = T(a 



and i/ien corresponding solution to (1.1) blows up infinite time in the regime described by Theorem 
lO 



The proof of Proposition 3.5 relies on a monotonicity argument on the energy £ which is the core 
of the analysis, see Proposition |4.6[ and the strictly outgoing behavior of the unstable mode induced 
by the non trivial eigenvalue — £ < of H , see Lemma 4.10 The fact that the regime described by 
the bootstrap bounds (3.31), (3.32), (3.33), (3.34), (3.35) corresponds to a finite blow up solution 
with a specific blow up speed will then follow from the modulation equations and the sharp derivation 
of the blow speed as in [3T] . 



4. Improved bounds 

This section is devoted to the derivation of the main dynamical properties of the flow in the 

The three main steps are first the derivation of a 



bootstrap regime described by Definition 3.4 



monotonicity property on £ which allows us to improve the bounds (3.31), (3.32), (3.33), (3.34) in 



[0, Ti(a+)], second the derivation of the dynamics of the eigenmode and the outgoing behavior of the 
unstable direction, and eventually the derivation of the sharp law for the parameter b which allows 



to bootstrap its smallness (3.31) and will eventually allow us to derive the sharp blow up speed. 



Remark 4.1. All along the proof, we will introduce various constants C(M), 5{M) > which do 
not depend on the bootstrap constant K{M). An important feature of all these constants is that, up 
to a smaller choice of b*{M) or a larger choice of K(M), we assume that any product of the form 
C(M) f(b) where lim^o fip) = or any ratio S(M)/K(M) is small in the trapped regime. This will 
be used implicitly in this section. 
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4.1. Coercitivity of £. Let us start with showing that the linearized energy £ yields a control of 
suitable weighted norms of (w,e) in the regime t G [0,Ti(a+)]. 

Lemma 4.2 (Coercitivity of £). There exists Mq > 1 such that for all M > Mo, i/iere exisi^] 
<5(M) > and C(M) < oo such that in the interval [0,Ti(a+)), there holds: 



'fwj 2 + 



(<9 r w 



,21 



C(M)[if(M)] 



i 6 4 

4 ■ 



|log6| 2 



(4.1) 



Proof of Lemma 14.21 

This is a consequence of the explicit distribution of the negative eigenvalues of H and the a priori 
bound on the unstable mode (3.35). Indeed, let t G [0, Ti(a+)), then first observe from (3.21), (3.22), 
d3T35b that 



|( e ,V0l 2 + |(<9 s e,v0l 2 < M 2 + M 2 + IM 2 (d 6 p Bl ,v>) 2 



< 



[tf(M)] 



i 6 4 

4 - 



+ C(M)6 2 |6 S | 2 < [K(M)]3- 



|log6| 2 ^ |logfe| 2 

where we used the estimates of Proposition |2 . 1 1 and the well localization of tp. This yields 



(4.2) 



^4(^,V'a) 2 + ^(^,Va) 2 



s [K ^m~ 2 



+6 2 



y 4 (l + |log(y)|) 2 



+ 



|Ve| 



(4.3) 



and similarly using the orthogonality condition (3.4): 
1 



-^(w, (xm$)\) 2 + ^(d t w, (xm$)a) 2 



= (6Ae,XM*) 2 
< b 2 M c 



y 4 (l + |log(y)|) 2 



+ 



|Ve| 



(4.4) 



Moreover, applying Lemma C.3 yields: 



A 



\He\ 



> 5{M) 



|Ve| 



+ 



y 2 y\l + \\og(y)\f 
Introducing the rescaled version (C.13) of Lemma C.3, we then conclude: 

|Ve| 2 



£ > 



1 



\\H x wY + 8x{M) 

2 



o 



£ 



y 4 (l + |log(y)|)< 



a / (vd tW y + 

Vel 2 ' 

+ " 



+ 



y4(l + |log(y)|)^ 



C{M)[K(M)} 



i b A 

4 - 



|log6|' 



> 



\ 2 {H x w) 2 + 5{M)\ 2 



{Vd t wf + 



(d r w) 



C{M)[K(M)\< 



|log6p 



where we used the Hardy bound (C.3), and fl4.1to is proved. This concludes the proof of Lemma 4.2 



Remark 4.3. Note that (4.1) together with the Hardy estimate (C.l), the coercitivity estimate 

iper: 

(4.5) 



(C.9) and (4.4) yield the following weighted bound on e which will be extensively used in the paper: 

iV-? + l 



let 



ri{s,y) = X 2 + d t w(t,\y) = d s e(s,y) + bAe(s,y), 



^recall remark 



4.1 
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then: 



y\l + |lo g2 /| 2 ) 



+ 



+ / |Vr?| 5 



< 



< 



c(M) 



\£\ + [K(M)]i 



iM)\E\ + y/K{M)-. 



|log6| 2 
b 4 



(4.6) 



(4.7) 



|logfo| 2 ' 

4.2. First bound on b s . We now derive a crude bound on b s which appears as an order one forcing 
term in the RHS of the equation for e (3.11 ). This bound is a simple consequence of the construction 
of the profile Qb and the choice of the orthogonality condition (3.4). 

Lemma 4.4 (Rough pointwise bound on b s ). There holds the rough pointwise bound^ 

{e,H$) 



b s + 



(4.8) 



(AQ,*) 

Remark 4.5. This is in contrast with |31] where the b s term could be treated as degenerate with 
respect to e thanks to a specific choice of orthogonality conditions and the factorization of the 
operator H in the wave map case. This difficulty in our case will be treated using a specific algebra 
generated by our choice of orthogonality condition (3.4) which gives the right sign to the leading 
order terms involving b s in the energy identity (4.6), see (4.24), (4.38). 

Proof of Lemma 14.41 

Let us recall th at th e equation for e in rescaled variables is given by ( |3.11 ), (3.12), (3.13). Observe 
also that from (1.19), the adjoint of H B with respect to the L 2 (R N ) inner product is given by: 

H* Bl = H Bl + 2b 2 D. (4.9) 
To compute b s we take the scalar product of (3.11) with xm&- Using the orthogonality relations 

(d?e,XM*) = {d?(P Bl - Q),XM$) = 0, V m > 
we integrate by parts to get the algebraic identity: 

b s [{AP Bl , xm$) + 2b(Ad b P Bl , xm<&) + (Ae, X M®)\ 

= -(* Bl ,XM$) - (£,H Bl ( X M$)) + 2b(d s e,A(xM$)) + (N(e), X M^). (4.10) 
We first derive from the estimates of Proposition |2.1[ 

6 4 



(^,XM^) 2 < |log6|2 . 



Similarly, using (4.6) yields: 

(d s e,A( X M^)) 2 <C(M) 

and 



c{M)\£\ + y/K(M) 



|log6p 



(4.11) 



(4.12) 



(e,H Bl ( X M<S>)) = (e, HQ) - (He, (1 - X m)$) + O M c b\ c(M)\£\ + y/K(M)- 



b 4 



logfop 



We then use the improved decay (2.2) and (4.7) to estimate: 

2 



(He,(l- XM m 2 < (I t^4v) 

\Jy>M 1 + V N J 



~ M 



b l 



logfep 



6 



recall remark 



4.1 
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Thus: 



b 4 



|log6|* 



(4.13) 



similarly, 

(AP Bl , xm$) + 2b(Ad b P Bl , XAf^) + (Ae, Xm$) 

= (AQ,$) + 

= (AQ,$) + 



log(6) 
b 

log(6) 



6 4 



(4.14) 



where we have used that in the trapped regime £ < K (M)b 4 / [log(b)] 2 . Finally, on the support of \M 
and for b < b^M) small enough, the term Q dominates in Q b = Q + b 2 T\. Hence, for the nonlinear 
term, we have from Sobolev and (4.7): 

=■2 ,.3 



\(N(e), X M$)\ 



< 



< 



+ 



[1 + 



C(M) 



£ + yjK{M)- 



log6| 2 



Injecting this together with ( |4.1l[ ), ( |4.12[ ), ( |4.13| ), ( |4.14[ ) into ( |4.10[ ) yields p^8|f|and concludes the 



proof of Lemma 4.4 



4.3. Global H 2 bound. We derive in the section a monotonicity statement for the energy £ which 
provides a global H 2 estimate for the solution. The monotonicity statement involves suitable re- 
pulsivity properties of the rescaled Hamiltonian H\ in the focusing regime under the orthogonality 
condition (3.11) and the a priori control of the unstable mode (3.35), which themselves rely on the 



positivity of an explicit quadratic form, see Lemma 4.7 



Proposition 4.6 (H 2 control of the radiation). In the trapped regime, there exists a function T 
satisfying 

tS 1^(1)^ (4.15) 



M ' v "^'llog&l 
and such that, for some < a < 1 close enough to 1, there holds: 



d f £ + F 



dt I A 2 (!- Q ) 



< 



A 3-2a 



y/K(M)- 



log&p 



(4.16) 



Proof of Proposition 4.6 
step 1 Energy identity. 
Let 



V(t,r) 



N + 2 

z___Lo N ~ 2 

N 



1 



W (r) = y 2 V{- x ), V(y) 



N + 2 
N-2 



Q N ~ 2 (y). 



We first have the following algebraic energy identity which follows by integrating by parts from 

'(d t w) 2 



1 

2cft 



V{d t wY + / (i? A «; 



+ wH\w 



+ / d t wH x F Bl . (4.17) 



recall remark 



4.1 
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We now use the w equation and integration by parts to compute: 



dtVwH\w 



d t Vw(F Bl - d tt w) 



d_ 
~dl 



dtVwdtw 



}"/ 



d t VwF Bl - / dtV(dtw) 



(4.18) 

duVwdtw (4.19) 



We next pick < a < 1 close enough to 1 and combine the above identities to get: 



1 d 

2X^di 



,2a 



(dtrW) 



V(8 t w) 2 + / (H x w) 2 - 2 / d t Vwdtw 



-Ry + R 2 + 



2ab 



dtVwdtw 



where R\ collects the quadratic terms: 



Ri 



ab 
T 
b_ 
A^ 



(d tr w) 



V{d t wf + / (H\w) 



+ 



d t V{d t wf 



A 2 



« f {dyV? -a fvr, 2 + a J (He) 2 + ^j\2V + y VV)r] 2 
J e(2V + y ■ W)AQ 



and i?2 collects the nonlinear higher order terms: 



R-2 



d t wH x F Bl - / d t Vw 



d tt Vwd t w (4.20) 
d t V(AQ) x w 

(4.21) 
(4.22) 



step 2 Derivation of the quadratic terms and treatment of the b s term. 
Let us now obtain a suitable lower bound for the quadratic term R\. The main enemy is the b s 
term which is order one in e and will be treated using a specific algebra generated by the choice of 
orthogonality condition (3.4). 

N 

Observe from H(AQ) = that (AQ/X) x (y) = (l/A)^(AQ)(y/A) satisfies: 

-A(AQ/A) A (y) - (l/X) 2 V(y/X)(AQ/X) x (y) = 0. 
Differentiating this relation at A = 1 yields: 

H$ = H(DAQ) = (2V + y VV)AQ. 
We inject this into the modulation equation (|4.8|) to get: 



(4.23) 



/( I & I b^ 1 
e(2V + y ■ VV)AQ = b 2 s {$, AQ) + |fe s |0 ( U + y/K(M) 



log6| ; 



We thus conclude using the sign 



($,AQ) > 



and (4.21), (4.8) that: 



Ri > 



A 3 



« J (dyT]) 2 + A(3 - a)V + • VF]?? 2 + a J (He) 2 + ci( 



+ o(M + y^M)_|_ 



(4.24) 



for some universal constant ci > independent of M. 
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step 3 Coercitivity of the quadratic form. 
We now claim the following coercitivity property of the quadratic form in r\ appearing in the RHS 



of (4.24) in the limit case a = 1, see Appendix [B 



Lemma 4.7 (Coercitivity of the quadratic form). There exists a universal constant cq > such that 
for all 7] S H^ ad , there holds: 



[dyrif + 



2V+-yVV 



V 2 >co f(d y r } ) 2 --[{n^) 2 + {r i ^) 2 } 



From a simple continuity argument, there exists < a* < 1 such that given < a* < a < 1, for 
all rj G H^ ad , there holds: 



(3-a)V+-yVV 



*7 2 >f /(M 2 -^[( 7? ^) 2 + (r ? ,$) 2 ]. 



a J (d yV ) 2 + 

We now pick once and forall such an a < 1 and control the negative directions 

b l 



Using (4.3) and (4.7), it yields: 



{T),ty 2 <b\S\ + y/K(Mj 



logfep 



Similarly, we compute (rj, <£) = (77, xm$) + (j], (1 — xm)^) for which (4.4) and (4.7) yield 

{V,XM^? <b\E\ + y/K[M)-^ 



and we have, applying (C.l): 

7 2 

>y>M/2 V 

This together with (4.24) yields the lower bound on quadratic terms 

_^ 
A3 



< ^ / 1^ 



i?i > 



6 4 



(4.25) 



c l {(b s ) 2 + \£\) + o(^mM) ][logh][2 
for some universal constant c\ > 0. Indeed, a straightforward integration by parts in (3.19) yields: 

£<j \dyV\ 2 + j \He\ 2 . 
step 4 Control of lower order quadratic terms. 



The lower order quadratic terms in (4.20) are controlled similarly 



dtVwdtw 



~ A 2 



1 + 



■ + 



< 



< 



h ( 6C(M)|f 1 + ^>s# 

A 2 " 



(s + ^%^>) • (4 - 26) 



and, with the help of (3.32), 

duVwdtw 



- \ A 3 A 3 



+ 



^[bC(M)\S\ + ^K(M) 



If 
b 4 



log6| 2 ) 
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Remark 4.8. We note here that (4.26) is sufficient for the proof of our theorem. Indeed, the es- 
timated term J dtVwdtw has been integrated by parts with respect to time, so that it becomes a 
part of J-. Furthermore, we note that to compute (4.16), we multiply T b y X 2a . Consequently, the 
commutator ba/X J dtVwdtw appears on the right-hand side. However, (4.26) yields that, in the 



og&| . Similar arguments 



trapped regime, this supplementary term is controlled by b/ \ 3 y/K(M)b 4 /\ 
will be repeated implicitly below for the terms which require an integration by parts with respect to 
time. 



step 5 Rewriting of the nonlinear R 2 terms. 



It remains to control the nonlinear R 2 terms in (4.20) given by (4.22). According to (3.17), this 
term contains b ss type of terms which cannot be estimated in absolute value and require a further 
integration by parts in time. Let 



F Bl =F X - d t F 2 with F 2 = j(d s P Bl )x 



(4.27) 



and rewrite: 
R2 



d t wH x F 1 



dtVw 



F 1 + ^(AQ) A 



d tW H x d t F 2 + / d t Vwd t F 2 



We now integrate by parts in time to treat the F 2 term: 



d t wH x d t F 2 



d t Vwd t F 2 

d t wH x F 2 - J d t VwF 2 



d 
dt 



{duVw + 2d t Vd t w)F 2 + J d tt wH x F 2 . 



The last term is rewritten using (3.17) and integration by parts: 



f d tt wH x F 2 = ([Fx- d t F 2 - H x w]H x F 2 



1 d 
2dt 



WFo 



VF4 



d t VFi + [F 1 - H x w]H x F 2 . 



eventually arrive at a manageable expression for R 2 : 



R-2 



d_ 
dt 



d+Vw 



d t wH x F 2 - 
F 1 + ^(AQ) X 



d t VwF 2 + 



1 



WF, 



VF4 



(4.28) 



+ / d t wH x F x - / (duVw + 2d t Vd t w)F 2 



We now aim at estimating all the terms in the RHS of (4.28). According to (3.17), we split Fx into 
four terms: 

Fx + ^|(AQ) A = [* Bl + Fx,x + F 1;2 + N(e)] x (4.29) 



A 2 



with 



Fx,x = bAd s P Bl + b s (AP Bl - AQ), Fx, 2 = [f'(Q) - f'(P Bl )] e. 



(4.30) 



step 6 Fx terms. 
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The F\ terms are the leading order terms. 

rom (2.9 

l*B!|< 



^> B x terms: We first extract from (2.9) the rough bound: 

b 2 



which yields: 



|log&|(l + y 2 
1 + |lo g2 /| 2 



+ C(M)6%< 2i?1 



(4.31) 



l*R I 2 < — 



and thus from (4.7): 



< 



< 



A 3 J (1 + y 4 ) 



b b 2 
A 3 |logfe| 



C(M)J\£\ + y/K(Mj 



b 4 



|logb| ! 



Next, we use the fundamental cancellation H(AQ) = and (|2.9[) to estimate: 



6 4 ri + |log(6y)l 1 , 

l2<y<2S + A 2„2|,-.. l i 1 «a. 



|log6| 
b 2 

+ (l + y *)|Iog6|^ fl1 / 2 ' 



&V|log&| ^<v<2Bi 



+ 



log(M) + |log(l + j/)| . 
l + y 2 



Ly<2B 1 



and thus 



Hence: 



(4.32) 



< 1 
~ A 3 



L 2 



6 2 



(l + y) 2 |ff(* Bi : 



Fx i terms: From (2.7), (2.8), there holds : 



l^i.il < 



l + |log(6y)| 1 , 



+ 



log(M) + |logy| 



|logbj *<v<^ 6 2 y 2 |log6| -f<y<^ 1 1 + y2 
and, recalling that differentiation w.r.t. y acts as a multiplication by 1/(1 + y) : 



\HFx,x\ <C{M)± 
from which 



\b s \b 2 



+ y 2 



l + |log(6y)| 



log(M) + | logy |" 



j(l + y 2 )\H(Fx,x)\ 2 <\b 



|lo g 6| l2 <s/<^ + 6 2 y 2 |log6| lB i<y<*Bi + \ + y 2 

(1 + |logy| 2 ) 



\\ogb\' 



-\Fi,x\ z <\b s \ 2 b 2 . 



(i + y 4 



(4.33) 
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Hence similar arguments as with the *&b 1 terms yield: 

b 



d t VwF 1 



b 4 



< ~b\b s \c(M)j\s\ + ^m llogb? 



and 



dtwHxFij 



< 



C(M)b \b s \ 
A 3 \\ogb\ 



\£\ + ^/K{M)- b ^< ° 



|log6| 2 ~ A 3 



|log6| |logfe| 



|log6|< 



F\i terms: The explicit expansion of the cubic nonlinearity and the bound (2.7) yield: 



l*1.2| < 



C(M)b 2 



i + y 2 



kl, |VFi 2 |< 



C{M)b 2 . C(M)b 



from which: 



1 

A 2 



d t Vw(F li2 ) 



C{M)b 3 



< 

~ A 3 



1 + y 3 
< b 



e + 



i + y 2 



We\ 



(4.34) 



1 + y 6 ~ A 3 

and, after integration by parts of the laplacian term: 



b\£\ + y/K(M) 



V 1 



\\ogb\< 



1 

A 2 



dtwH x (F 1>2 ) 



< 



C{M) 



b 2 , b 2 , , 

\e\ + rrlV^I 



+ r 



< A 

~ A 3 



Nonlinear term N(s): We expand the nonlinearity: 

iV(e) = 3P Bl e 2 + e 3 . 



This yields using (3.27), (C.l) the rough bound: 



£ 



\N(e)\ < , ; - 



In what follows, we will use the following bound on r\ which follows from (4.6), (C.l): 

b 4 N 3 



IOTI 



L-<|V7?| L2 < c(M)|£| + V^(M)- 



log6| 2 



We then estimate: 
1 



A 2 



9tV«;(JV(e))A 



< A 

~ A 3 



/ T+7 ~ a7 |V£|l2 ( c(m)|£:| + 



log&p 



for 6q < b*(M) small enough. We split the second term: 



d t wH x 



A 2 



Vd t w ■ V 



A 2 



Vdtw 



(N(e)) : 
A 2 



(4.35) 
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The second term is estimated in brute force: 



Vdtw 



(N(e)h 
A 2 



< 



1 



< 



A 3 J 1 + y 5 ~ A 3 
1 



\vv\ 



i + y 6 

3 



< 



b 6 4 



A 3 |log6| 2 ' 



The first term in (4.35) is split into two parts: 



/ 



Vd t w ■ V 



V8 t w • [V(w 3 ) + 3(P Bl ) x V(w 2 )] + ^ / e 2 Vr/ ■ VP Lh . 



The last term is integrated by parts in space and then estimated in brute force: 



A 3 



s 2 V V ■ VP Bl 



3 

A 3 
1 



< 

~ A 3 



J v [e 2 AP Bl + 2eVP Bl ■ Ve 
e\\Ve\ 



1 + y 4 + 1 + y J 



~ A 3 



+ 



1 



The first term is the most delicate one and requires first a time integration by parts: 
j Vd t w ■ [V(w 3 ) + 3(P Bl ) x V(w 2 )] = j t {J \Vw\ 2 ^w 2 + 3(P Bl ) x w 
- 3 J wd t w\Vw\ 2 -3 j \Vw\ 2 [wd t (P Bl )x + (P Bl )xd t w] . 



We may now estimate all terms in brute force: 
"3 



\Vw\ 



-w z + 3(P Bl ) x w 



1 f \Ve\ 2 1 6 4 

< -g[mL«> + \yPB 1 \L^\\ys\Lo a J < 



wdtw\VwY 



< -T2\ye\L°°\yr)\L° 



< 



b b A 

A 3 |log&| 2 ' 



|Ve| 



2 ^(< M )l f l + ^Wn^2 



/ 



\Vw\ 2 wd t (P Bl ) : 



< 



\yw\ L oo f \Vw 



A 3 



y 



l + y< 



< -^\Ve\ L ^l + C(M)\b s \ — 



■+C(M)b\b S \ly< Bl 

\\ogb\\ f\Ve\ 2 



< 



b b 4 
A 3 |log6|< 
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where we used the rough bound extracted from (2.8): |%PbJ < C(M)bl y < Bl , and finally: 



\Vw\\P Bl ) x d t w 



< 



< 



1 

b b 4 

A3 |log6| 2 : 



iv^i 2 < 

i + y 3 ~ 



C(M)£ + y/K{M) 



b 4 V 



|log6| s 



for 6o < b*(M) small enough. The above chain of estimates together with remark 4.8 closes the 
control of the nonlinear term N(e). 

step 9 F2 terms. 



We estimate from (2.8): 



and hence: 



dbPBl + /' V ^' 2 ^' /tt^ Pbi 

l^ll^l + J \v v \\vd b p Bl 



(i + y) 

d t wH x F 2 



~ |log6| 2 



(4.36) 



< JM 

~ A 2 



(1 + y 4 



< 



1 

A 2 |log6| 



3(M)|5| + 



b 4 



\logb\< 



< 1 
~ A 2 



dtVwF2 



< 



\b s \b f \d b P Bl \\e\ < |6 fl |6 



A 2 J (l + y 4 ) ~ A 2 |log6| 



:(M)|£| + v 7 K(M) 



|log6| 2 



< 1 
~ A 2 



/ |VF 2 | 2 + 

similarly: 

/" (^tVta + 2d t Vd t w)F 2 + 



< 



A 2 



\dbP Bl \ 
(l+y 4 



+ J \vd b p Bl \ 



d t VFl 



< JM 

~ A 3 

< 1 
~ A 3 



< 1 (6 a ) 2 < 1 b 4 
~ A 2 |log5| 2 ~ A 2 |log6| 2 ' 

f((|6 s | + 6 2 )|e| + %|)|d 6 P Bl 

(N + &) 



< J6- 
~ A 3 



(i + y 4 



+ \b s \b 



\d b P Bl 

1 + y 4 
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|log&| 



(M)\£\ + yZpZ)-^ + 



|log&| 2 |log&| 2 ' 



Eventually, (4.32), (4.33) ensure: 



b 6 b 2 \b s \ 



|log6| 2 |log6| 2 



< 

~ |log6| 2 
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which together with (4.36) yields: 



< 1 b 3 \b s \ < b 



A 3 |log6| 2 ~ A 3 llogbl 1 



We similarly estimate from (4.34) and after integration by parts: 

b 2 \e\\8 b P Bl \ 



I 



{Fx,2)\H\F2 



A 3 



J l + + y |Wflil viTF + i 



6 2 |Vd 



< 



C(M) 



A 3 |log6| 



1 + 



+ 



+ r 



|Ve| 2 \ 5 < 6 6 4 



A 3 |log6| 2 ' 



For the nonlinear term, we extract from (|2.8l) the rough bound 



\H(dbP Bl )\ < [C(M)+log(b)]—^l y < Bl , 



which together with (C.l) ensures: 



1 

Y 2 



(N(e)) x H x F 2 



< [C(M)+log(6)] 



A 3 

6 s ||log6| 4 



Z — : n — l?y<Bi 

l + y 2 l + y y " 



£ A3 



(l+y 4 )|logy| 2 



< 



A 3 



>/6(c(M)|f| + \/K(M)-. 



6 1 



< 



log6| 2 J ~ A 3 |log6| 2 ' 



step 10 The remaining i 7 ^ term has the right sign. 
It remains to estimate the term 

f H x wH x F 2 



in the RHS of (4.28). Let us stress onto the fact that this term is a priori no better O(-^S) due to 



the b s contribution and the bound (4.8), recall remark 4.5 



We now claim that the main contribution has the right sign again. 



Indeed, we first compute from the T\ equation (2.16): 



HT\ = -$ + c bX BoAQ, HdtTi = O 

4 



1 1 2<^<^ 

6|tog6| (1 + y 2 ) 



(4.37) 



We then apply the decomposition (2.22): 



H(d b P Bl ) = h(wTi + 2b{ XBl ~ + b 2 d h \og{B 1 ) PBl T l + b 2 XBl d b Ti) = -26$ + E 



and estimate using (2.8), (2.21), (4.37): 



i + y 2 



1 



|lo g 6| ^<v<^ + b 2 y 2 \logb\ "tP<» 
I) 2 



In particular, 



|log6| 
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and thus using the modulation equation ( |4.8[ ) : 

A 3 



H x wH x F 2 



2^b s (e,H<S>) + ^0 



(He)H(d b P Bl ) 
1 6,1 



A 3 



J He (-26$ + E) 



6 1 



|log6| 2 



'A 3 



(AQ,$) 



+ o 



( £ ,F<I>) + -0 



|log6| s 



26 (e,F$) 5 
M 



£\ 
M 



+ A 3 ° V |lo«6| 2 



(4.38) 
(4.39) 



The recollection of all above estimates yields (4.16) and concludes the proof of Proposition 4.6 



4.4. Improved bound. We now claim that the a priori bound on the unstable direction (3.35) 
coupled with the monotonicity property of Proposition 4.6 imply the following improved bounds: 

Lemma 4.9 (Improved bounds under the a priori control ( |3.35[ )). There holds in [0,Ti(a+)]: 

b(t) 



\\(Vw(t),d t w(t) + jt^{(1 - XBx(b(t)))^Q))x(t)\\i^xL^ ^ 6 |log6 |, 
b 4 (t) ^ 6 4 (0) 



> 



|log6(t)| 2 A 2 ( 1 - Q )(t) " |log6(0)| 2 A 2 ( 1 -«)(0) ' 
K(M) 6 4 



|6,| 2 < 

W)\ < 



2 |log6| 2 ' 

K(M) 6 4 
2 (log6) 2 ' 



(4.40) 
(4.41) 
(4.42) 
(4.43) 



Proof of Lemma 14.91 

step 1 Energy bound. 



The energy bound (4.40) is a consequence of the conservation of the energy. Indeed, the conser- 



vation of the energy and the initial bounds of Lemma |3.1| ensure 

E{u, d t u) = E(u Q , ui) = E(Q) + O(6 0V /|log6 | 



(see Appendix 3.1) and thus: 

£(Q) + O(6o|log& |) 
= \J [dt(P Bl )x + d t w] 2 + \ J \V(P Bl ) x + V w ?- l -j [{P Bl ) x + 



{AAA) 



w 



We lower bound the first term by expanding 
d t (P Bl )x + d t w 



d t w + h -{{\ - xbJAQ)a + 1(xbAQ)x + jiHXB.mx + j(d b P Bl )x 



d t w + j((l-x Bl )AQ) x + E, 
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with 



£ 2 <bg|log6o|, 



where we used the bootstrap bounds (3.31), (3.32). Finally : 

b 



I 



A 



((l- XBi )AQ) x + d t w 



O(b 2 \logb 



[dt{P Bl )x + d tW ] 2 > 1 - J 
We then expand the second term: 

\ J [V(P Si )a + VH 2 ~\f \{P Bl )x + wf = \ J [VP Sl + Ve] 2 -\J [P Bl + ef 



(4.45) 



|VP Bl | 



(4P Sl e 3 + e 4 ) 



\P Bi \*-(e,AP Bl +p- Bi ) + 



We\ 



3 / Pie 2 



From the construction of Pg 1 , 

\ J \VP Bl \ 2 -\ J \P Bl \ 4 = E(Q) + 0(b 2 \logb\). 



(4.46) 



The linear term is treated using (2.9), the improved decay (2.2) and (4.31): 
\(e,AP Bl + P Bl )\ = 



\(e,b z DAP Bl -V Bl )\ 

< \\e/y\\ L 2\\y(b 2 DAP Bl - * Bl )\\v < b\Ve\ L 2 

We now rewrite the quadratic term as a small deformation of H and use the coercivity bound ( C 
to ensure: 



(4.47) 



J \Ve\ 2 - 3 J P 2 Bl e 2 > co J \Ve\ 2 + Def, 



(4.48) 



with 



Def := 3 



J{Q 2 -p 2 B y 



co 



Collecting ( [2.7] ) and (C.l), on the one hand, and ( |4.2| ) on the other hand, we compute: 
j (Q 2 - P 2 Bl )e 2 < \\y 2 (Q 2 - P 2 ^ ||V E || 2 2 < 6||V E || 2 2! \(e,^)\ 2 < 6 2 |logfe|. 
The nonlinear term is easily estimated from Sobolev: 

J \(3P Bl +e)e 3 \ < \\yP Bl \M\ye\\L~\\Ve\\ 2 L2 < \/b~ \\V e\\ 2 L2 
Injecting ( ggg ), ( pL47] ), ( g46| , ( [449] ), ( [448] ), ( [450] ) into gig yields flOo] ). 
step 2 Lower bound on b. 



(4.49) 



(4.50) 



We now turn to the proof of (4.41). First observe from the bootstrap estimate (3.32) that 



b 2 



1 — a , 



This implies: 



d 



46 3 



(4.51) 



1 — a . 
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and ( |4.41| follows. 

step 3 Improved H 2 bound. 

We now turn to the proof of ( |4.43 ). We integrate (4.16) in time and conclude from (4.1), (4.15): 



\£{t)\ 



VA(o). 

+ {K{M))l 



b\t) 



(4.52) 



|log6(i)|- 



+ [A(t)] s 



(!_«, /•* b(r) b\r) 
J [A(r)]3-2a |l og 6(r)P 



We then derive from ( 4.51| ): 

"* b(r) 6 4 (r) 



< 



< 



o [A(r)]3-2«|l og 6(r)|2 
1 b\t) 



dr 



4 

A 3-2a | log6 |2 



2(1 -a) A 2 ( 1 -«)(*)|log6(i)| 2 2(1 -a) Jo A 3 ~ 2a |log&| 2 
* 4 (*) , HTT^f b(r) b\r) 



|log6| s 



A 2 ( 1 -«)(i)|log6(t)| 2 
and hence the bound: 



+ y/K(M) 



[A(t)] 3 - 2q |lo g 6(r)| 2 |log6(r; 



/o [A(r)]3-2» |logt(r)|2 ~ |logK()P' 

Injecting this into (4.52) and using the initial bound ( A.12 ),( A.17| ) and the monotonicity (4.41) 
yields: 



S(t) 



< 



|logfe(0)P 



|log&(t)| s 



< y/K(M)- 



b\t) 



logb(t)|5 



(4.53) 



and (4.43) follows. (4.42) now follows from \AA) and (4.53) 



This concludes the proof of Lemma 4.9 



4.5. Dynamic of the unstable mode. We now focus onto the dynamic of the unstable mode. We 
recall the decomposition 



Y(t) 



(9 s e, if>) 



(4.54) 



and the variables given by (3.22): 



«+(*) = M s ) + ^^(dbPBi,^), «_(a) = 2_(s) - ^= 



Lemma 4.10 (Control of the unstable mode). There holds: for all t £ [0, Ti(a+)], 



and k + is strictly outgoing: 



a7c_l_ 



log&| 



(4.55) 
(4.56) 
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Proof of Lemma 14.101 

We compute the equation satisfied by the unstable direction (e, ip) by taking the inner product of 



(3.11) with the well localized direction ip to get: 

d 2 



ds 2 



- CM) = E{s) - {d A s P Bl ^) 



(4.57) 



with 



E(e) 



-- -(* Bl ,ip) -b s (AP Bl ,i>) -b(d s P Bl +2Ad a P Bl ,i/>) -b(d s e + 2Ad s e,^) 

b s (Ae,iP) + (N(e),iP) + b 2 (Ae,Dip) + ((f'(P Bl )-f'(Q))e,iP). (4.58) 

Simple algbebraic manipulations using (4.54), ( |3.22 ) and the initial condition yield the equivalent 
system: 

d r~ . E + (s) d r~ . , E-(s) , . 

with 



(4.59) 



E+(s) = E(s) - |(5 6 P Bl ,^) £L(«) = tf(a) + ^(d b P Bl ,^) 



(4.60) 



We now have from the explicit formula (4.58), (4.60), the exponential localization of ip, the orthog- 
onality 

(V,AQ) = 0, 



the estimates of Proposition 2.1 and the bootstrap estimate (3.32) the bound: 

b 2 



\?\E±\ < \b\(\b s \ + ^\ + ^K(M)-t-) < Vb 



|log6| ' 



|logfe| 



(4.61) 



which together with (4.59) yields (4.56). Let then 



Q = K 2 



|logfr| 2 



then from (4.59), (4.61), (3.32), we estimate: 

ck_ llogfel 2 9 , 

jiog6p r / r 



ds 



4|log6| 2 21og& 



fe 5 



+ 



VC|log6| 2 2 |log6| 2 /- 6 2 VC r ^. 



We integrate this in time 



Q{s) <g(0)e-^ s + I e~^ s -^da<\ 



where we used the initial inequality (A. 18) yielding that £7(0) < 1. This concludes the proof of (4.55) 
and of Lemma 14.101 
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4.6. Derivation of the sharp law for b. We now turn to the derivation of the sharp law for b 
which will yield the required monotonicity statement on b to close the smallness bootstrap estimate 



(3.31), and will eventually lead to the derivation of the sharp blow up speed (1.10). 



Lemma 4.11 (Sharp derivation of the b law). Let 

Pb = x^oQi 



G{b) = b\AP Bo \ 2 L2 + / b(d b P Bo ,AP Bo )db 



Z(«) = (9 s e, AP Bo ) + b(e + 2Ae, AP Bo ) + b s (d b P Bo , AP Bo ) - b s (^(P Bl - P Bo ), AP Bo 
then there holds: 

G(b) = 646|log6| + 0(b), |Z| < K(M)b, 



d_ 

ds 



{G(6)+Z(s)} + 326 s 



K(M) 



(4.62) 

(4.63) 
(4.64) 

(4.65) 
(4.66) 



Remark 4.12. Observe that (4.65), (4.66) essentially yield a pointwise differential equation 

b 2 



2|log6| 

which will allow us to derive the sharp scaling law via the relationship 
Proof of Lemma 14.111 



b. 



The proof is inspired by the one in |31j . We multiply (3.11) with AP Bo and compute: 



b s AP Bl + b{d s P Bl + 2Ad s P Bl ) + d 2 s P Bl ,AP Bo ) = -(* Bl ,AP Bo ) - (H Bl e,AP Bo ) 
8 2 s e + b(d s e + 2Ad s e) + b s Ae, AP Bo ) + (N{e), AP Bo ) 

We further rewrite this as follows: 

(b s AP Bo + b(d s P Bo + 2Ad s P Bo ) + d 2 P Bo ,AP Bo ) = -($ Bl ,AP Bo ) 

- (b s A(P Bl -P Bo ) + b(d s (P Bl - P Bo ) + 2Ad s (P Bl - P Bo )) + 8 2 (P Bl - P Bq ),AP Bo ) 

- (H Bl e, AP Bo ) - (d 2 s e + b(d s e + 2Ad s e) + b s Ae, AP Bo ) + (N(e), AP Bo ). (4.67) 
We now estimate all terms in the above identity. 

step 1 b terms. 



An integration by parts in time allows us to rewrite the left-hand side of (4.67) as follows: 

d r 



,AP Bo +b{d s P Bo +2Ad s P Bo ) +d 2 P Bo , AP Bo ; fj 



G(b) + b s (d b P Bo ,AP Bo ) +\b s \ 2 \d b P Bo \ 2 L2 (4.68) 



with G given by (4.63). Observe from ( 3.32[ ) the bound 



IM 2 |^PB li 2 <^<(^(M)) 2 -^< 



b 2 



b\ 2 - ^gb\' 
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We now turn to the key step in the derivation of the sharp b law which corresponds to the following 
outgoing flux computation^ 



(* Bi ,APb ) = 326 2 1 + 



|log&| 



as 6—^0. 



(4.69) 



Indeed, we first estimate from (2.9): 



(*Bi -c b b 2 X BoAQ,AP Bo ] 



< b q 



l + |log(6y)| , l + |log(l + y)| 



+ 



|log6|(l + y 2 ) (l + y 2 ) 2 



< 

~ |logb| 



The remainder term is computed from (|2.10|) and the explicit formula for Q (|1.3|): 

1 



{ebb Xb^AQ, AP Bq ) 



1 + 



2|log6| 
32b 2 ( 1 + O 



|log6| 
1 

log&| 



{KQY + 0(l) 



and (4.69) follows. 

We now estimate the lower order terms in b which correspond to the second line of (4.67). One term 
is reintegrated by parts in time: 

-(d 2 (P Bl - Pb ),ap Bo ) = ~ [b s (d b (p Bl - p Bo ),ap Bo )} + bUd^PB, - p Bo ),d b AP Bo ). 

The remaining terms are estimated in brute force using (2.8) and (3.32) which yield: 



s A(P Sl - P Bo ) + b(d s (P Bl -P Bo ) + 2Ad s (P Bl - P Bo ))),AP Bo ] 



+ b 2 s 
step 2 e terms . 



{d b {p Bl -P Bo ),d b AP Bo \ 



\b I 2 b 2 



b 2 



log&| 



We are left with estimating the third line on the RHS of (4.67). We first treat the linear term 



from (4.1), (4.7), (3.34): 



/ 



(H Bl e,AP Bo ) < \(He,AP Bo )\ + / \e\\P z Bl - Q 2 \\AP Bo \ + b 2 (DAe,AP Bo ) (4.70) 



On the one hand, (4.7) together with bootstrap estimates yield: 



\e\\P 2 Bl -Q 2 \\AP Bo \ < b 2 



; / j — 



y<B 



(i + y 2 ) 2 



:i + y) £ 



< JL 

~ |log(* 



see again |31| for more details about the flux computation statemement and its connection to the Pohozaev inte- 
gration by parts formula 
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On the other hand, after integration by parts, we repeat the same arguments and apply (C.4). This 
yield: 



(DAe,AP Bo ] 



< V 



< &i 



6 7 2 

+ 6 2 



y<B 



£ \ ,2 

+ b 2 



(1 + V 4 ) JB /4<y<B /2 (1 + V 2 ) Jy<B 1 + V 



|Ve| 



y 



(1+1/5) 



+ 



+ 



B /A<y<B /2 (1 + y A ) J \Jy<B 1 + V 



|Ve| 



V^loi(fe)| (c(M) |£ | + s/K(M) 



< 



vl lo g( 6 )l 



llogftl 1 



Finally: 



|(ife,AP Bo )| < |tfe| L2 ^|Ioi&7+0^M)- 



logbp 



We further integrate by parts in time to obtain: 
(d 2 e + b(d s e + 2Ad s s) + b s Ae, AP Bo 

with 



d_ r 
ds 

6 S 



(d s e,AP Bo ) + b(e + 2Ae,AP Bo 
(d s e + bAe,Ad b P Bo ) + (e, $ ft 



= -AP Bo - A 2 P Bo - bAd b P Bo - bA 2 d b P Bo . 



We thus estimate from (|4.lD, (|4.5|), (|4.7|), (|3.32|), (|3.34|): 

16. 



(5,e + 6Ae, Ad b P Bo ) + (e, $&) 



< 16, 



+ 



?<y<Bo 2/ Jv<Bo 1 + 2/ 2 



6 4 6 2 

^~^ (M) ^ii67 



The non linear term is estimated as previously. Indeed, we have: 
(N(e),AP Bo ) < /(|P Bl | + |e|) e 2 |AP Bo | 



< ^l|y(|i'B 1 | + N)||L<-||(l + I/ 2 )APBolU« 



B ° lei 2 



y(i + y 4 ) 



< ^(M) 
6 2 



5 + JT(M)- 



log6| 2 



K(M) 



step 5 Control of G(6) and X. 



Injecting the estimates of step 1 and step 2 into (4.67) yields (4.66). It remains to prove (4.65). 
The estimate for G(b) is a straightforward consequence of the choice (4.62) and the explicit formula 
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Ql -3p . It remains to control X. We integrate by parts in space in (4.64) to rewrite: 

Z{s) = (d s e + bAe, AP Bo ) - b(s, AP Bo + A 2 P Bo ) + b s (d b P Bo , AP Bo ) - b s (d b {P Bl - P Bo ),AP Bo ) j . 
The b terms are estimated as in step 1: 



|6 S | (d b P Bo ,AP Bo ) - (d b (P Bl - P Bo ),AP Bo ] 



< i-^i < b. 
~ b ~ 



The linear term is estimated using (4.1), (4.5), (|4.7h, (3.32), (3.34): 



(d s e + bAe, AP Bo ) - b(e, AP B() + A 2 P Bo 



< 



+ b 



< I 
~ b 



12X1 



+ 



6 2 



y<B y 4 ( 1 + l lo §yl 



K(M)b 



< 



and (4.65) is proved. 



This concludes the proof of Lemma 4.11 



5. Sharp description of the singularity formation 

We are now in position to conclude the proof of Proposition |3.5| and Theorem |1.1| as a simple 
consequence of the a priori bounds obtained in the previous section. The proof relies on a topologi- 
cal argument which closes the bootstrap argument, and then the sharp description of the blow up 
dynamic is a consequence of the a prori bounds obtained on the solution and in particular the mod- 



ulation equation (4.66). 



Proof of Proposition 3.5 



We argue by contradiction and assume that for all a + G 

Ti{a + ) < T{a + ). 



|log&o| ' |log&o| 



(5.1) 



Lemma 4.10, a simple continuity argument ensures that Ti(a+) is attained at the first time t where 

\^(t)\ ~ 

The fundamental fact now is the outgoing behaviour ( 4.56| ) which together with (5.1) ensures 



dt 



> 0. 



Thus from standard argumentj^J the map 



3_ 



|logbo| ' | logbo I 

a. 



H- T 1 (a + ) 



is continuous. 



We may thus consider the continuous map: 

$ : 



|log&o|' llogboM 
a + 



K+{l 1 {a+)) 6 2 (Tl(o+)) 



'see Lemma 6] for a complete exposition 



34 



M. HILLAIRET AND P. RAPHAEL 



On the one hand, (5.1) implies: 



h 2 



|logh | ' |logfe | 



C {-1,1}. 



On the other hand, the outgoing behavior (4.56) together with the initialization k+(0) = a+ ensures: 



h 2 



|logbo| 



-1, $ 



h 2 



log6 | 



1 



and a contradiction follows j^jThis concludes the proof of Proposition 
Proof of Theorem 11.11 

step 1 Finite time blow up and derivation of the blow up speed. 



3.5 



Let from Proposition 3.5 an initial data with Ti(a+) = T(a+). We first claim that u blows up in 
finite time 

T = T(a + )<+oo. (5.2) 

Indeed, from (4.41), 

A 2(l-a) < 6 3 d th A | < x 2 -^ < b = _ A 

Integrating this differential inequation yields 

*<As(0)-A3(t)<l 



and ([5^ follows. The (H l C\H 2 ) x (tfnH 1 ) bounds ( ]3.33[ ), ( |3.34| ) on (e, S^e) and hence on (it, 5 t «) in 
the bootstrap regime and standard H 2 local well posedness theory ensure that blow up corresponds 
to 

A(t) ^0 as i-> T(a+). 

We now derive the blow up speed by reintegrating the ODE (4.66) and briefly sketch the proof which 
follows as in [3T] . 

First recall the standard scaling lower bound 

A(t) <C(u )(T-t) 
which implies that the rescaled time is global: 

"* dr 



s(t) 



Let 



so that from (4.65 ): 



l 7 = 646|log6| 1 + 



1 



|log6| 



o A(r) 

J = G + 1 

ie b 



+oo as t — > T. 



J 



64|logJ| 



and J satisfies from (4.66) the ODE: 

Js + 



J' 1 



128|lo gl 7| 2 



l + O 



|logJ| 



1 + 



0. 



(5.3) 



'This topological argument is of course the one dimensional version of Brouwer's fixed point argument used in [5]. 
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lloe Tl 2 

We multiply the above by 1 j$ 1 , integrate in time and obtain to leading order: 

1 \ \ . A, 21ogs / „ / 1 



IgQggg / +0 



ie _^ =i , = 21 ^ l + ( — 
Vllogsl// A s V V v/|logs| 



where we used (5.3). Integrating this once more in time yields: 

1 



-logA = (logs)" 1 + O 



\/ 1 logs | 



and thus 



b = -A t = exp -V|logA| 1 + O 



|logA| 



Integrating this from t to T where A(T) = yields the asymptotic 
A(t) = (T — t)exp ( -V|logA(t)| ( 1 + O 



1 



|logA(t)| 



which yields (1.10) 



step 2 Energy quantization. 



It remains to prove (1.9) which can be derived exactly as in [31], this is left to the reader. 
This concludes the proof of Theorem |1.1[ 



Appendix A. Modulation theory 



This appendix is devoted to the proof of Lemmas 3.1 and |3.3| The arguments are standard in the 
framework of modulation theory and we briefly sketch the main computations. 



A.l. Proof of Lemma 13.11 First note that the bounds 

|V(P Bl - Q)\v + b\AP Bl - 6(1 - XBjAQltf < 6|log6| 

ensure that our initial data are of the form 

uo = Q + fjo, ui = fji 

for a small excess of energy in the sense that: 

\Nvo,m\\ L 2 xL 2 <b \logb \, ||V 2 ?7o, Vr?i|| L 2 xL 2 < 6 , (A.l) 

Hence the continuity of the flow associated to (1.1) ensures the existence of a time To > (uniform 
in fjo,iji) for which the solution u to (1.1) (uo,ui) satisfies on [0, To]: 

sup \\V(u-Q),d t u\\ L 2 xL 2 < fo |logfoo], - (A. 2) 

10,1b] 

Step 1 : Modulation near Q. 

The non degeneracy (AQ,&) ^ ensures^] that u admits on [0, To] a decomposition 

u(t) = (Q + e(t))x(t) (A.3) 

with: 

(e(t),XM^) = 0. (A.4) 



ii 



as a direct consequence of the implicit function theorem and the smoothness of the flow (1.1 1 
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Moreover, A G C 2 ([0, To]; M+) and noting that % satisfies 



we obtain the bound: 



l A <°> - !l s nth- 

|log£>o| 

We then let 6(t) = -At(t) on [0,T ]. 

Step 2 : Positivity of 6. 
Straightforward computations yield: 

d t i(t)=(d t u-^Au) . 

Taking the scalar product with ^m$, we obtain at the initial time: 

((ui) i ,XM&) 



(A.5) 



6(0) = A(0) 



MB) 



((A«o) i ,XM$)' 

A(0) 



(A.6) 



where (2.5) together with (A.5) imply: 



(W_L,XM*) = 6 (AQ,xm$) + O . 

MO) Vl lo g(M 

((Auo) i ,XM*) = (AQ,XAf*) + 0(^|log(6o)|). 

A(0) 



(A.7) 
(A.8) 



This yields the positivity of 6(0) and moreover, the positivity of b(t) for small time together with: 



h 2 



|log(6 ) 



(A.9) 



b(t) = b + O 

As 6 > 0, we may introduce the decomposition: 

u(t) = (Q + e) x(t) = (P Bl {b(t)) + e)\(t) ie e(t) = e(t) - (^(6(0) " Q)- ( A -!0) 
Observe from §L4§ , ( |A.4[ ) that 

Vie[0,T o ], (e(i),XAf*) = 0. (A.ll) 
The uniqueness of such a decomposition is guaranteed by the (local) uniqueness of (A,e). 

Step 3 : Smallness of e. 

To complete the proof, we obtain the smallness of e in H 1 and H 2 . To this end, we note that: 



e(0) = (uo)_i_-P Bl (b(o)) 



P 



Bx(b )Ji ~ P Si(6(0)) 

A(0) 



A(0) 



Simple computations based on the estimates of Proposition |2.1| yield the expected result : 



iVe 



I? 



< 



6o|log(6 )| 



e(0) 



i + y 4 



+ llV 2 e 



L 2 ^ 



L 2 



log(6 )| 



(A.12) 
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A. 2. Proof of Lemma 3.3, The proof of this lemma is divided into two steps. First, given 
(rjQ,r]i,d + ) satisfying smallness condition (3.1) for small 6q, we prove that b,b s and w satisfy (3.31)- 
(3.34). Then, we show that, given (60,770,^1), we can apply the inverse mapping theorem to 
d+ i — y k+(0) close to 0. The arguments are standard and we refer to [5] for a detailed proof in 
a similar setting. 

Step 1: Smallness of initial modulation given (770, 771, 



Given (770, 771, <Z+) satisfying smallness condition (3.1) we can apply Lemma 3.1 this yields To and 

(A.13) 



b,e,w such that (3.31) holds and 

llVio 



L 2 ^ &o|log(&o)| 



\V 2 w 



68 



log( 



We emphasize in particular that Lemma 3.1 implies bo/2 < 6(0) < 2bo for sufficiently small bo- 

As previously, we focus now on bo unds satisfied initially. We first compute b s (0) using (1.1) 
and the orthogonality condition (A. 11). Recalling that (d^P Bl , X M$) = = f° r an Y 

integer k, we get like for (4.10): 

b s [(AP Bl , xm®) + 2b(Ad b P Bl , XM$) + (Ae, X M$)] 

= -(¥ Bl ,XM*) - (^^(xm^)) + 6(5.e,A(xM*)) + {N(e), X M$) 

where, denoting LHS and RHS the left-hand and right-hand side at initial time, we compute, for 
sufficiently small bo w.r.t. M : 



\RHS\ < C(M) 



h 2 



log( 



+ \\ d s e \\L 2 (y<M) ) ; 



16,(0)1 



(AQ,xm$) < \LHS\. 



On the other hand, after time-differentiation, we obtain : 



d s e(0) = A(0)$e(0) = -b s (0)d b P Bl{m) - b(0)Au + A(0) (b AP Bl{bo) ) , 

A(0) 



(A.14) 
(A.15) 



Observe now from (2.8) that 



l^ P Bi(6o)l 



L 2 {y<2M) 



< C(M)b < Vbo 



which together with (A.5), (A. 9) and (3.1) yields: 



l^(0)|| L 2 (y<2M) ~ A(o)||a t£ (o)|| L 2 (y<2M) ^ |7— 7j 



+ \b s (0)\y%, 



(A.16) 



which together with (A.14) concludes the proof of the initial bound (3.26) on b s . 

6(0) 



Then, we compute : 



d t w(0) = ui 



V A(0T l(6(0)) Aid) /;i 



A(0) 



so that, introducing (A.15) and previous estimates on 6(0), we compute : 

®(a-YR,™™UQWnJr.a < 6 n | ln(6 n )l < y/bo 



11^(0) + ^((1 - Xb i( 6(0)))AQ)a(q)II^ < 60I ln(6 )| 



and 



\\Vd t w(0)\\ 



L 2 < 



65 



log6 | ' 



(A.17) 



Together with (A.13), this yields (3.27) and (3.28). 
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Finally, straightforward computations yield: 

1 



^(d s e,i>) - ^ (d b P Bl ,ip) 



Consequently, we apply (3.28), noting that w(t) = (e(t))x(t), and (A. 15) because of the exponential 
decay of tp to compute 



M0)I<i4t 

|log£»o| 



(A.18) 



Step 2: Computation of d+. 

We now claim from an explicit compuation that given a+, the initialization (3.24) can be reformulated 
in the form 



dF 

F(d + ) = a + with ^-|d + =o 



2 

L- 



O(b 



(A.19) 



dd + ,a+ ~ w 2 

which from the implicit function theorem concludes the proof of Lemma |3.3| 
Let us briefly justify (A.19). We want to study the mappping 

d + .— ► [b(t),b s (t),(e(0),i>),(d s e(t),ip)] 

where V is a neighborhood of 0. To this end, it is necessary to study the dependencies of all initial 
parameters on d+. For conciserness, we denote by d differentiation w.r.t. <i+ in what follows 
Computation of (A(0), e(0)). As a first step in the modulation theory, we proved that (A(0), e(0)) = 
3>(«o) where $ is a smooth mapping H X (R N ) -> M x F 1 ^) defined in a neighborhood of Q. Due 
to the exponential decay of tp £ C°°(M. N ) we thus have that A(0) is a smooth function of d+ with 
differential d\(0) = dX G K. We have the same result for e with differential de(0) = de £ i/ 1 (M iV ). 
By definition, we have 

e(0) = n - Qi 

so that: 



A(0) 



i . 

A(0) 



Computation of 6(0): From (A.6), 6(0) is a C 1 mapping with: 

i2 



db(0) 



r/A 



A(0) 



+ 



((«i) i 

-HQ) 

((Au )^,xa/$) 

A(0) 

((«i) i ,xm$)((AV>) i ,XM$) 

MO) A(0) 

((An ) i ,xm$) 2 

A(0) 



((A 2 u )^,Xa/$) - ((A«i) i ,XAf^) 

A(0) A(0) 



((A«o) i ,Xa/$) 2 

A(0) 



where (A.6) and (A. 7) ensure that, for some db £ M, there holds : 

db(0) = db + O(b ). 

Computation of e(0): Next, 

e(0) = e(0) - (P Bl (6(o)) - Q) 
Consequently, (e(0),ip) is also a smooth function of d + with derivative dpsi(0) satisfying 

dpai(0) = (de,if>) - d6(0)(^P Bl(6(0)) , V) 

Replacing de by its values, and applying that (AQ, tp) = together with |A(0) — 1| < 6§/ |l°s(^o) I ? we 
get: 

(de.VOHMI^ + OC&o) 
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so that: 



d pSl (0) = U\\ 2 L2 +0(bo). 



Computation of d s e(0) + b s (0)dbPB 1 (b(o)) : From (A. 15), 
d s e(0) = -b s (0)d b P Bl{m) 



6(0)An + A(0) (b AP Bl(bo)) 



A(0) 



so that (5 s e(0) + b s (Q)d b P Bl ( b r \\, V 1 ) is a smooth function of d + with derivative : 



dps 2 (0) = -db(0)(Au ,il>) + dX 



(b AP Bl(bo) ) + {b A 2 P Bl{bo) )_^_ 



4 -6(0)(A^), 



where, for the same orthogonality reason (AQ,^) = 0, we have: 

(An ,^) = (AQ,^) + 0(6o) = 0(&o) 

Consequently (fps2(0) = 0(6q)- 



Conclusion: Finally, there holds 

1 

2 



K+(0) 



1 



(e(0),V) + ^=(3 s e(0) + 6.(0)^P fll(6(0))s ^) 



and k+(0) = a + reduces to a simple ID equation F(d+) = a + with F computed as combination of 
the above functions so that it is smooth in a neighborhood of 0. Moreover, there holds: 



dp*i(0) + -^=dps 2 {0) 



+ O(b ), 



and (A. 19) is proved. This concludes the proof of Lemma 3.3 



Appendix B. Coercivity estimates 

The aim of this section is a proof of the coercivity properties of the quadratic form: 

B(rj,rj) = (Bv,v) = [ \d r r]\ 2 + f Wrj 2 

where 

3 fi Q r 2 

W (r) = 2V+^rV 



2 (l + r 2 /8) 2 4(l+r 2 /8) 3 



is a 

2 



We use the elementary method developed in [8]. The coercitivity property of Lemma 4.7 
consequence of the two following facts. First the index of B on H} = {u radial with J |Vn| 2 + J \ < 
+00} is at most 2. From standard Sturm Liouville oscillation theorems, see Theorem XIII. 8 [33], 
this is equivalent to counting the number of zeroes of 

f BU = on (0,oo), 

(B.l) 

1 U(0) = 1 U'(0) = 0, 

and this can be analytically reduced to counting the number of zeroes of a Bessel function. Then 
we need to show that the orthogonality conditions (77, ip) = (77, $) = are enough to treat the 
two negative directions. Arguing exactly as in [S], see also [13], this is equivalent to first invert 
the operator B on H} ad , and then show that B restricted to Span{B~ 1 ip,B~ 1 &} is definite nega- 
tive, which is an elementary numerical check. We shall check these two facts below and refer to 
[8] for the proofs that this implies the claimed coercitivity property. Note that the proofs in [8J are 
given for exponentiallly decaying functions and potentials, but one checks easily that the decay of the 
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1- 
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it) 40 60 80 100 120 140 160 180 200 



-1- 



Figure 1. Solution to (B.l) computed by MAPLE 



potential |W(r)| ~ \ at infinity and |3>(r)| ~ \ are more than enough to have all proofs go through. 



B.l. Computation of the index of B. We claim: 

Lemma B.l (Derivation of the index). The index of B on is at most 2.. 
Proof. First, we note that W{r) > W{r) where: 

W(r) 



2(l + r 2 /8) 3 ' 

Hence, classical Sturm-Liouville theory ensures that U has less zeros than U the unique solution to : 



r 3 fa 



dr 



+ WU = on(0,oo), 



(B.2) 



U(0) = 1 U'(0) = 0, 
Second, we look for U of the form: 

U(r) = ^U(r 2 /2), 

with U a sufficiently smooth function. Denoting by s the new variable r 2 /2, straightforward calcu- 
lations yield that U is a solution to : 

' d 2 - 

^U + WU = on(0,oo), 

ds 2 V ' h (B.3) 

I U(0) = U'(0) = 1, 
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where: 



Setting then U{s) = y/l + a/4 U(l/y/l + s/4), we obtain that J7 is a solution to (|B.3l) if and only 
if U is a solution to 

( d 2 ~ d ~ 

r 2 ^C/ + r— C/+ (96r 2 - 1)17 = on (0,1), 
dr dr 

[ [7(1) = l>(l) = -8, 

Hence, U is a combination of Bessel functions: 

U(t) = CiJ(l,4>/6r) + C 2 F(l,4V6r) 

We compute (Ci,C2) and draw the explicit combination with MAPLE. We obtain Figure [2j The 
computed solution U has two zeros on (0, 1). Moreover, it diverges in so that U(r) ~ K/r close to 
with K 7^ As a consequence 

if 

C/(r) ~ — when r — > oo, 

and thus the index of — A + W on H^ ad is exactly two. Hence the index of B is at most 2. This 
completes the proof of Lemma B.l □ 



0,25- 



0,2- 



0,15- 



0,1- 



0,05- 



0,0- 
-0,05- 
-0,1- 



I | I | I | IM I | I I^K | I | I | I | 
0,1 0,2 0,3 0,4^—6^—^6 0,7 0,8 0,9 1,0 
tau 



Figure 2. Solution to (B.3) computed by MAPLE 
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B.2. Choice for the orthogonality conditions. We now invert B. We first check numerically 
that the solution U does not vanish at infinity ie 

lim U(r) > 0, 

r— >+oo 

see Figure [TJ 

Hence U is not a resonance -note that if U had been a resonance, we could have removed the 
resonance by diminishing a bit the potential and getting a potential with index 2 and no resonance-, 
and thus from standard ODE arguments, [3], there exists unique smooth solution in H l ra d of: 

' <! " ' " +WU = ip on (0,oo), 



BU 



r 3 dr 
U'(0) = 0, 



4u" 

dr 



with (l + r^){7 G L c 



(B.4) 



and 



fit/ 



d 



^ 77 
dr 



+ W?7 = $ on (0,oo), r 2 

V ; with 1 + )U e L c 

logr 



(B.5) 



C/'(0) = 0, 

We denote B~ 1: ijj and B~ l § the respective solutions to these systems. We recall the explicit formula 

In the remainder of this section we, check numerically that the restriction of B to Span(/3 -1, 0, B~ l Q) 
is definite negative, or equivalently: 

Lemma B.2 (Numerical check of the orthogonality conditions). The symmetric matrix 

' (B^ifj,^) (B- 1 ®,^) ~ 



satisfies: 

and is thus definite-negative. 



(rTV,VO<0 and detB>0, 



(B.6) 



Numerical proof of Lemma B.2 We use standard MATLAB routines for the computation of 



solutions to ( B.4|B.5 ). We note that we only fixed the initial value for U'(0). The value U(0) is left 
open in order to achieve the expected decay at infinity which characterizes the inverse. In order to 
obtain B~ l ijj, we first compute ip. We obtain that the corresponding eigenvalue is approximatively I = 
—0.5860808922. Because ip decays exponentially, we only need to obtain an approximation on a short 
time-range. We computed our solutions until T^^ max = 30. We emphasize here that we use an explicit 
scheme. As a drawback, the accumulation of errors tends to make the numerical solution to become 
negative when the exact solution is exponential small. Hence, our scheme becomes unstable after 
time T^^ max = 18. Nevertheless, we extend our numerical solution with after this time. This induces 
an exponentially small error. The pictures in Figure [3] illustrate this computation. On the left-hand 
side is drawn the obtained solution. On the right-hand side, we draw iptest{f) = '0(r)exp(\/— lr). We 
observe here that our solution enters the exponential asymptotic regime before the instability comes 
into play. 

The solution B~ l ip is computed with the extension of ip. Straightforward ode analysis shows that 
the unique solution decaying fast at infinity behaves like 1/r 2 asymptotically. The choice of U(0) is 
made with respect to this criterion. Figure [4] illustrates that we obtained a solution with the suitable 
decay. As previously, on the left-hand side is a picture of the numerical solution. On the right-hand 
side we plot B~ l iptest{ r ) = r 2 i3 _1 V'(r). In the latter computations, this solution is involved in scalar 
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products with ip. Hence even if drawn until T 7 
solution until Tg-i^ max = 18. 



APsi on r=0..3O 

0.2 I . . 1 




-0 1 
-3 2 
-0.3 

-0.4 

~° 5 50 1 00 1 50 200 250 300 



= 300, we only need a precise computation of this 



Asymptotic behavior of Apsi 
50 | 1 1 1 • 




50 100 150 200 250 



Figure 4. Numerical simulations for B ip 

The last solution B^ 1 ^ is computed with the same method. In this second case, the expected 
decay of the solution is log(r)/r 2 . Figure [5] illustrates that we obtained a solution with the suitable 
decay. The picture on the right-hand side restricts to the time-interval r = 0..100 because this is the 
significant region. In the latter computations, this solution is involved in integrals which converge 
slowly. Hence, we compute this solution until max = 1000. 

We now compute numerically the entries of the matrix IB . We first compute (B~ 1 Q,?p) = {B~ l ijj,&). 
The exponential decay of tp implies that we need to compute the first integral ip) on a shorter 

time-interval. Hence, we prefer this computation to the second one. We compute the L 2 -scalar 
products with a standard trapezoidal method. Changing the time-interval and the time-step, the 
computations are stable up to an error of 10~ 2 . We get the following approximations for the integrals 
involving ip. 

(B~V, VO = -4.63 ± 10" 2 (B' 1 ®, V>) = 32.65 ± 10~ 2 . 
The last integral is a more involved computation. Indeed, standard real analysis implies that there 
holds : 

I(M) := / i3 _1 $Q(r)$(r)r 3 dr = (B^ 1 ®, $) + err(M) 
Jo 
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Asymptotic behavior of AO 

100 | 1 1 1 




Figure 5. Numerical simulations for B 



with a remainder satisfying err(M) = (K + o(l)) ln(M)/M 2 for some constant K. This remainder 
going slowly to 0, we see numerically that our computations has not converged even when integrating 
until T B ~i^ <max = 1000 (see Figure |6| red crosses). In order to improve the rate of convergence we 
compute an approximation of coefficient K and substract the estimated error term of our computa- 
tions. This yields Figure [6j blue circles. On this second computation we obtain a very good rate of 
convergence. Hence, we provide the approximation 

(B" 1 ®, $) = -574.25 ± 10~ 2 

Hence 

det(B) = 1591 ± 10 



which concludes the numerical proof of Lemma B.2 



Appendix C. Some linear estimates 

We start by recalling some obvious integration-by-part results : 
Lemma C.l. For any N > 3, there exists a constant C for which there holds, for any v 6 d (M. N ) 

\v(y)\ 2 ~ 1 



+ sup \y\ 2 \ V (y)\ < C 

Irn \y\" 

Looking for control on further derivatives, we prove: 



\Vv(y)\< 



Lemma C.2 (Hardy inequalities). Let N = 4. Then \/R > 2, Vi> G H? ad (R N ) 
following controls: 



IVP < 



y < R y4(l + |l og y|)2 ~ y v < B y 



:a,) 2 , 



+ 



R<y<2R U 



<logi? 



1^1 

y<R V 2 



+ 



V<2 



y<2 



(C.l) 



there holds the 



(C.2) 
(C.3) 

(C.4) 
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similarly, using f(y) 



j/ 3 (l-log(j/)) 

\v\ 2 



e<y< 



a y 4 (l - logy) 2 
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we get: 

M 2 v-/ 



a<y<R 



< 



1 - log(y) 
v(a)\ 2 + 



+ 2 



vd,,v 



y<a 



y 3 (l - logy) 

i 

2 



y < R y 4 (l + \\ogy\ 



\dyV\ 



y<R V 



(C.7) 



(C.5), ( C.6| ) and (C.7) now yield (C.3). The last inequality (C.4) is a straightforward variant of [5T] 
Lemma B.l, (B.4)] and is left to the reader. □ 

Lemma C.3 (Coercitivity estimates with H). Let ip be the first eigenvector of H . Then there exists 
c > and Mq > 1 such that for M > Mq, there exists 5{M) > such that given u £ H^ ad (R N ), 
there holds 

(Hu, u)>cf (d y uf - - [(«, V) 2 + (u, XM®?] (C.8) 



/ 



(/7-u) 2 > 8{M) 



(d y uf 



y 



+ 



2/4(1 + |logy|) 2 



1 



5(M) 



(U, XM^f 



(C.9) 



Proof. (C.8) is a standard consequence of the coercitivity of the linearized energy which admits 
exactly ip as bound state and AQ as resonance at the origin, the good enough localization of $ (2.1 ) 
and the nondegeneracy (2.2). The detailed proof is left to the reader. 
To prove (C.9), we first observe the key subcoercivity property: 



J (Hu) 2 = J (An + Vu) 2 = j (An) 2 -2 J V(d y u) 2 + J (AV + V 



> c 

where we used the asymptotic value 

N(N + 2)(N -2) 



i + y 4 



1+y* 



(CIO) 



V{y) 



IT 



i + o(^) 



as y —¥ +oo. 



(C.9) now follows by contradiction. Let M > fixed and consider a sequence u n such that 

,2 



"7 



y 4 (l + |logy|) 2 



1 



and 



[{Hu n ) 2 <-, (u n ,XM$) = 0. 



(C.ll) 



(C.12) 



Then by semicontinuity of the norm, u n weakly converges on a subsequence to G solution 
to Huqq = 0. n^ is smooth away from the origin and hence the explicit integration of the ODE and 
the regularity assumption at the origin Uqo E H} oc implies 

Uoo = aAQ. 



On the one hand, the uniform bound (C.ll) together with the local compactness of Sobolev embed- 
dings ensure up to a subsequence: 



(9yU n ) 

i + y 4 



+ 



\U-r, 



1 + 



(dyUpp) 

i + y 4 



+ 



1 +2T 
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thanks to the xm localization. We thus conclude that 

a(AQ, Xm&) = (uoo, Xm$) = and thus a = 0. 



On the other hand, the subcoercivity property (C.10), the Hardy control (C.2), (C.3) and (C.ll), 

(d y u n ) 2 



(C.12) ensure 



i + y 4 



+ 



"7 



> C > 



from which 



a 



{d y KQf 

i + y 4 



+ 



|AQ| 2 
1 + 



(dyU C 



i + y 4 



> C > and thus a / 0. 



A contradiction follows. This concludes the proof of (C.9) and Lemma C.3 



□ 



Straightforward computations show that the coercitivity estimates with H can be adapted to any 
of the operator H\ yielding, for any A > and u G H^ ad (M N ), 



(H x u, u)>cj {dyuf - -L [(«, (^)a) 2 + («, (xm$)a) 2 ] 
for the same c and <5(M) as in Lemma C.3| 



(C.13) 
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